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Abstract 

It was found recently that processes of multidimensional tunneling are generally 
described at high energies by unstable semiclassical trajectories. We study two obser- 
vational signatures related to the instability of trajectories. First, we find an additional 
power-law dependence of the tunneling probability on the semiclassical parameter as 
compared to the standard case of potential tunneling. The second signature is substan- 
tial widening of the probability distribution over final-state quantum numbers. These 
effects are studied using modified semiclassical technique which incorporates stabiliza- 
tion of the tunneling trajectories. The technique is derived from first principles. We 
obtain expressions for the inclusive and exclusive tunneling probabilities in the case of 
unstable semiclassical trajectories. We also investigate the "phase transition" between 
the cases of stable and unstable trajectories across certain "critical" value of energy. 
Finally, we derive the relation between the semiclassical probabilities of tunneling from 
the low-lying and highly excited initial states. This puts on firm ground a conjecture 
made previously in the semiclassical description of collision-induced tunneling in field 
theory. 
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1 Introduction 



Tunneling in systems with several degrees of freedom is an exceptionally rich subject of 
investigation IH [2]. The features and probability of multidimensional tunneling depend 
crucially on the properties of underlying system, or rather on the degree of complexity of 
its classical dynamics. In particular, expressions for the tunnel splittings of energy levels 
are qualitatively different in the cases of integrable [21 IH IS] and near-integrable [SI El El IH] 
dynamics. The other drastically different case, tunneling in irregular (chaotic or mixed) 
systems, has been a subject of continuous theoretical [TOl [TTl [121 [IS El ESI HSl HZ] and 
experimental [HI [HI [20l [21] research for the last few decades. 

The basic concept in multidimensional tunneling is dynamical tunneling [23] . It is 
related to the classical dynamics and reflects the fact that transitions of a multidimensional 
system between the in- and out- regions of phase space may be classically forbidden even if 
there is no energy barrier separating the regions. In this case the quantum probability V of 
transition is on general grounds exponentially suppressed, 

V = Ae-^/'\ (1) 

where F and A are the suppression exponent and prefactor respectively. The transition 
itself is called dynamical tunneling [23], since the reasons for its exponential suppression are 
hidden in the particularities of classical dynamics. 

A new mechanism of dynamical tunneling has been independently discovered in Refs. [211 
[25] and [26] . It governs tunneling in non-separable systems with multiple degrees of freedom 
at energies exceeding certain critical energy Ec- The value of the latter energy depends on 
the details of the system dynamics but is always greater than the height of the potential 
barrier between the in- and out- states of the process. The new mechanism is general: it is 
relevant for tunneling in regular f2E{ [271 [28j and irregular [211 US] scattering problems, for 
transitions in time-dependent one-dimensional potentials [231 [29l [30] , in the case of chaotic 
tunneling [T2l [31] . Another example emerges in field theory where the new mechanism is 
generically inherent in the processes of collision-induced tunneling at high energies [32] [33]. 

The defining characteristics of the new mechanism have been given within the semiclassi- 
cal approach. It was noted that the semiclassical trajectories describing tunneling transitions 
acquire qualitatively new properties at E > Ec- Instead of connecting directly the in- and 
out- regions of the process, the trajectories end up performing unstable motion on the bound- 
ary between the regions. In the simplest case of two degrees of freedom this unstable motion 
proceeds along the periodic orbit describing oscillations on top of the saddle point of the 
potential. Following Ref . [2^] , we call the latter orbit sphalerori^ (or simply unstable periodic 
orbit) . 

^In chaotic case the new mechanism imphes anomalously weak falloff of particle wave function in some 
parts of classically forbidden region ("plateau structure" [l2l[24j). This behavior leads to anomalously large 
tunneling probabilities, the effect known as chaos-assisted tunneling [lOj . 

^This term is standard in field theory [34]; it is based on classic Greek adjective aLpaXepo^ — "ready to 
faU." 
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In general case of systems with more than two degrees of freedom the boundary between 
the in- and out- regions is normally hyperbolic invariant manifold (NHIM) p!S]; the trajec- 
tories in the new mechanism get attracted to this manifold. In this case we use the term 
sphaleron in the sense equivalent to NHIM. 

Due to the above property of the trajectories, tunneling at E > Ec proceeds in two 
stages. First, the long-living sphaleron "state" gets created. Second, the sphaleron decays 
into the final asymptotic region with the probability of order one. The overall transition 
remains exponentially suppressed, since creation of the sphaleron costs exponentially small 
probability factor. We call the overall transition sphaleron- driven tunneling. 

The aim of the present paper is twofold. First, we analyze the possibility of direct 
experimental observation of the the mechanism of sphaleron-driven tunneling. To the best 
of our knowledge, such observation has not been performed so far. We study two signatures of 
the new mechanism which may be helpful in future experiments. Second, we systematically 
develop modified semiclassical method for the calculation of tunneling probability in the 
sphaleron-driven case. 

We discuss two experimental signatures of the new tunneling mechanism. In Ref. |28] we 
have found that the probability of sphaleron-driven tunneling contains additional power-law 
dependence on h as compared to the ordinary case of potential tunneling. The additional 
factor is h}^'^ in the case of inclusive tunneling processes, i.e. processes without specification 
of the out-state. In this paper we review the result of Ref. p8] and extend the analysis to 
the new case of exclusive processes, i.e. processes with fixed out-state quantum numbers. 
We show that the additional factor is h in this case. For example, consider two-dimensional 
inclusive processes. Then, the prefactor A in Eq. ([T]) is proportional to h}^"^ and h in the 
cases of potential and sphaleron-driven tunneling respectively. For exclusive processes this 
dependence is h {h^) in the potential (sphaleron-driven) case. 

It is important to stress that the dependence of the tunneling probability on h can, in 
principle, be studied experimentally. Indeed, the semiclassical parameter, which we denote 
by h for convenience, is in fact a certain dimensionless combination of the Planck constant 
and parameters characterizing the system. Changing the latter parameters one varies the 
value of effective h. 

The second manifestation of the new mechanism is spreading of the out-state of the 
tunneling process over an anomalously wide range of quantum numbers. This effect was 
originally observed in Ref. [25] in the case of a one-dimensional system with time-dependent 
potential; here we show that it is present in the multidimensional case, cf. Ref. ^D]. Phys- 
ically, the widening of the out-state is related to the fact that the intermediate sphaleron 
orbit is classically unstable; thus, classical trajectories describing sphaleron decay spread 
exponentially over phase space. In quantum case this corresponds to the final state wave 
function which is almost constant in some region of quantum numbers. 

In the second part of this paper we develop the modified semiclassical technique which 
is essential in the case of sphaleron-driven tunneling. The motivation for the new technique 
becomes clear if we try to apply the standard method of complex trajectories to the prob- 
lem of inclusive sphaleron-driven tunneling in the scattering setup. Since the overall time 
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interval of the scattering problem is infinite, one generically finds two different trajectories 
corresponding to the two stages of the new tunneling mechanism: one trajectory starts in 
the in-region and tends to the sphaleron orbit as t — > +00, and the second trajectory starts 
at the sphaleron at t — *■ —00 and arrives into the out-region. The first of these trajectories is 
unstable: it can be destroyed by infinitesimally small changes in the initial Cauchy datal^l It 
is problematic to find unstable trajectories numerically. Besides, one wonders how to join the 
two trajectories in order to describe the overall two-stage process. Finally, it is not clear how 
to calculate the prefactor A of the tunneling probability. Indeed, the standard formula for 
the prefactor deals with the linear perturbations above the tunneling trajectory. When the 
trajectory in question is unstable these perturbations grow exponentially with time. Then 
the standard formula gives A = 0, which is obviously incorrect. 

Our modified semiclassical method overcomes the above difficulties. The main idea of 
the method was proposed in Refs. [2Sl ES]; here we present its detailed derivation. The 
modified method is summarized as follows. We evaluate the Feynman path integral for the 
tunneling amplitude in two steps. First, we restrict the integral to paths which arrive into 
the out-region in a fixed time interval r. Second, we integrate over r. The integration at 
step 1 can be done by the standard saddle-point method, since all trajectories at finite r are 
stable and interpolate between the in- and out- regions. On the other hand, the ordinary 
integral over r at step 2 should be evaluated with care. In particular, we find that in the 
case of sphaleron-driven tunneling this integral is saturated in the region r — > +00, rather 
than at the saddle point at finite r. 

The above manipulations with the path integral lead to a notably simple semiclassical 
description of sphaleron-driven tunneling. Namely, we show that the constraint in the path 
integral leads to the deformation of the semiclassical equations of motion with the imaginary 
term proportional to the small parameter e = e(r). The evaluation of the integral over r 
corresponds to taking the limit e —>■ +0 in both cases of stable and unstable trajectories. 
However, the resulting expressions for the tunneling probability are different in the two cases, 
since the integral over r is saturated in two different regions. In particular, the probability 
formula in the case of sphaleron-driven tunneling involves additional factor ^^/^ mentioned 
above. We call the modified semiclassical technique by the method of e-regularization. 

The new mechanism of tunneling is relevant only at sufficiently high energies, E > E^. 
Below Ec transitions proceed via the ordinary mechanism of potential tunneling. In ac- 
cordance with our results, the semiclassical expression for the prefactor A changes discon- 
tinuously across the critical energy. In particular, in the inclusive case in two dimensions 
A oc h}^'^ and h aX E < Ec and E > E^ respectively. This implies that both expressions break 
down in a small vicinity of Ec-, where the correct uniform approximation should be invoked. 
In the present paper we derive the required formula, which is continuous and applicable 
in the entire energy range. At \E — Ec\ 3> h^^"^ this formula coincides with the respective 
"potential" and "sphaleron-driven" semiclassical expressions. In this regard it is similar to 
the uniform approximation [9] for the tunnel level splitting at the point of transition from 
integrable to near-integrable systems. 

^Below we always refer to this kind of instabilities. 
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Next, we study semiclassically exclusive tunneling processes in the splialeron-driven case. 
We find that the new mechanism leads to proliferation of complex trajectories describing a 
given exclusive process. These trajectories form an infinite sequence and have the following 
structure: they get attracted to the sphaleron orbit, follow it for an integer number of 
periods and then slide away. The tunneling amplitude is the sum of the contributions of 
all these trajectories. In analogy to the case of inclusive probability the sum is saturated 
by the trajectories which spend an infinite time at the sphaleron. It is worth noting that, 
in contrast to the inclusive case, the individual trajectories describing exclusive process are 
stable. Thus, a priori, there is no need for the modified semiclassical technique in the case of 
exclusive transitions. Still, in this paper we demonstrate that our modified technique turns 
out to be useful in finding and organizing the tunneling trajectories. It also provides the link 
between the semiclassical descriptions of inclusive and exclusive processes. 

Finally, for the sake of completeness we study the processes of tunneling from low-lying 
in-states. Naively, such states and hence the corresponding tunneling processes cannot be 
described semiclassically; still, we show that the probabilities of these processes are given 
by the semiclassical formula ([1]). In addition, we show that the suppression exponent and 
prefactor of tunneling from the low-lying states can be obtained as certain limits of the 
corresponding quantities in the case of highly excited states. The limiting relation for the 
suppression exponent is known in field theory as the Rubakov-Son-Tinyakov conjecture [36]; 
it plays an important role in the semiclassical description of collision-induced tunneling [37j . 
We prove this conjecture in quantum mechanical setup. Our limiting formula for the prefac- 
tor shows that the probability of tunneling from the low-lying states contains a factor 
as compared to the case of highly excited in-states. 

We illustrate our findings by considering tunneling transitions in a simple model with 
two degrees of freedom. For this model we compare predictions of the modified semiclassical 
technique with the exact quantum mechanical results. The latter are extracted from the nu- 
merical solution of the stationary Schrodinger equation. We find perfect agreement between 
the two sets of results. 

The outline of the paper is as follows. After presenting the model in Sec. [2] we summarize 
the experimental signatures of sphaleron-driven tunneling in Sec. [31 In Sec. IHwe introduce 
the modified semiclassical technique: we review the standard semiclassical method in Sec. 14. 
introduce e-regularization in Sec. I4.2l and derive the uniform formula in Sec. 14. 3[ Application 
of the modified semiclassical method to the exclusive tunneling processes is discussed in 
Sec. O Finally, we study tunneling from low-lying in-states in Sec. El Section [7] contains 
discussion. Technical details are described in appendices. 



2 The model 

We start by introducing the scattering model of Refs. [3H1 ES] • It will be used throughout 
the paper for illustrative purposes. The model describes motion of a particle with unit mass 
in the potential 

y(a;,y) =cuV/2 + e-("+J')'/l (2) 
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Figure 1: (a) The contour plot of the potential (dashed lines) and the real part of the 
tunneling trajectory ed. E = 1.3, Ey = 0.05 (solid line). The saddle point is marked by 
the thick dot. (b) Time evolution of Rex for two complex trajectories with Ey = 0.05 and 
different values of total energy, E = 0.9 and 1.3. Note that Ec{Ey = 0.05) ~ 1.1. 

The potential represents two-dimensional harmonic waveguide extended along the x direction 
and intersected at an angle by the potential barrier. The contour plot of the potential is 
shown in Fig. [1^. In this and other figures we use the value ut = 1/2 for the waveguide 
frequency. Note that potentials similar to ([2]) typically arise in the studies of collinear 
chemical reactions [22] . 

We are interested in tunneling transitions of quantum particle between the asymptotic 
regions x — oo and x +oo of the potential (in- and out- regions respectively). In the 
in-region the particle evolves with constant momentum in the x direction oscillating along 
the y axis. The corresponding in-state is fixed by the total energy E and the energy of 
y oscillations Ey. Similarly, the out-state can be fully characterized by E and E^, where 
E^ is the final oscillator energy. In what follows we will often omit the specification of the 
out-state and consider the total (inclusive) probability of tunneling into the region x — >• +oo. 

The height of the potential barrier separating the in- and out- regions is Vq = 1. It is 
given by the value of the potential at the saddle point (x, y) = (0, 0). At < Vq the classical 
transitions between the regions are forbidden energetically, and their underlying mechanism 
is potential tunneling. On the other hand, it is shown in Ref. |38| that classical over-barrier 
transitions between the asymptotic regions take place at E > Eh{Ey), where Eh{Ey) is larger 
than Vq. Hence, at intermediate energies Vq < E < Ei,{Ey) the transitions are in the regime 
of dynamical tunneling, which we are interested in. 

As we have already discussed in the Introduction, the multidimensional processes of 
dynamical tunneling, such as ours, generically proceed via sphaleron-driven mechanism at 
sufficiently high energies. Let us illustrate the new mechanism in the model ([2]) comparing 
the behavior of semiclassical solutions at low and high energies [2S]. Consider the inclusive 
tunneling transition from the state \E,Ey) into the out-region x — >■ -|-oo. We postpone the 
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consistent formulation of the semiclassical method till Sec. HI The only fact we need here 
is that any tunnehng process is specified by a certain complex trajectory — solution to the 
(complexified) classical equations of motion. The latter should interpolate between the in- 
and out- regions of the process. 

Fig. [T]d shows the complex trajectories describing tunneling transitions at Ey = 0.05 and 
two values of total energy, E = 0.9 and 1.3. [The real part of the trajectory with E = 1.3 is 
also depicted in Fig. [T^.] The behavior of the two trajectories is drastically different. While 
the low-energy trajectory interpolates between the asymptotic regions x — *• ±oo, the solution 
with E = 1.3 gets stuck at finite x approaching the unstable periodic orbit as t — > +oo. The 
latter orbit is precisely the sphaleron discussed in the Introduction; it describes oscillations 
around the saddle point of the potential, see Fig. [1^. Clearly, the high-energy trajectory 
of Fig. [1] describes only half of the transition process, since it does not arrive into the 
out-region. Trajectory corresponding to the other half can be obtained by adding to the 
unstable periodic orbit infinitesimally small momentum in the direction of the out-region 
and evolving the system classically. Thus constructed, the overall semiclassical evolution 
involves two trajectories which describe creation and subsequent decay of the sphaleror0. 
This evolution corresponds to the mechanism of sphaleron-driven tunneling. 

One finds [2n] that there exists the critical value E = Ec{Ey) of total energy which 
separates the regions of qualitatively different behavior of tunneling trajectories. Namely, 
the trajectories interpolate between the in- and out- regions at E < Ec{Ey) and approach 
the sphaleron orbit at Ec{Ey) < E < Eb{Ey). This means that the mechanism of transition 
changes from potential to sphaleron-driven tunneling as the energy crosses the critical value. 
From the physical viewpoint Ec{Ey) can be understood as the energy of "phase transition" 
between the two regimes of tunneling. We remark that the energies of sphaleron orbits and 
thus the critical energy for the sphaleron-driven tunneling exceed the height of the potential 
barrier. Therefore, the new mechanism is relevant only in the case of dynamical tunneling. 

The region Ec{Ey) < E < Ei,{Ey) corresponding to the sphaleron-driven tunneling in 
the model ([2]) is shown in Fig. [2J The value of E(.{Ey) is found numerically by computing 
the complex trajectories at different energies and investigating their stability. 



3 Experimental signatures 

In this Section we show that the mechanism of sphaleron-driven tunneling leads to two 
observable effects which in principle can be used for identification of the new mechanism 
in future experiments. Both effects are related to the fact that the relevant semiclassical 
solutions are unstable. We illustrate our findings in the model ([2]) using the exact quantum 
mechanical results. The exact calculations of this and the subsequent sections are based on 
the numerical solution of time-independent Schrodinger equation, see Refs. [Ml [ISl EH] for 

"'There is another way to visualize the semiclassical evolution [25]. One introduces stable and unstable 
manifolds of the sphaleron orbit. These are formed respectively by the trajectories arriving at the sphaleron 
&i t ^ +00 and trajectories starting from it at i ^ — oo. Then, the evolution describing sphaleron-driven 
tunneling is guided in turn by trajectories belonging to the stable and unstable manifolds of the sphaleron. 
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Figure 2: Regions in the plane of in-state quantum numbers corresponding to the potential 
and sphaleron-driven tunneling mechanisms. 

the numerical method and Fortran 90 code. 

The first signature of sphaleron-driven tunneling is the direct consequence of the semi- 
classical analysis which will be presented in Sec. |H We find that the sphaleron-driven 
mechanism changes the power-law dependence of the transmission probability on h com- 
pared to the case of potential tunnehng. To be concrete, let us discuss inclusive tunneling 
transitions in the model ([2]). Then, the prefactor A of the probability is proportional to h)-/"^ 
and h in the cases of potential and sphaleron-driven tunneling respectively. 

The physics behind the additional power-law suppression becomes clear if one uses the 
qualitative analogy with the classically allowed creation of unstable state. The latter process 
considered at the classical level requires fine tuning of the Cauchy data. As a consequence, 
only a small part of the in-state wave function contributes into the amplitude of the process. 
This results in the additional suppression of the probability. On general grounds one expects 
similar formal suppression in the case of sphaleron-driven tunneling. 

Experimentally, one can try to observe the unusual power-law dependence on h by ana- 
lyzing the probability graph V{h). Note that the value of the semiclassical parameter which 
we denote by h is, in principle, adjustable in experiments. Indeed, the magnitude of quan- 
tum fluctuations is measured by the dimensionless ratio of the Planck constant to a certain 
combination of parameters characterizing the system. Changing the latter parameters in an 
appropriate way, one alters the value of the semiclassical parameter h without affecting the 
classical dynamics of the system. 

To illustrate this point consider the system ([2]). The key quantity which enters into 
the semiclassical expansion is the ratio of the action of the system to the Planck constant. 
Restoring the dimensionful units we obtain 



where ho stands for the physical Planck constant. In terms of dimensionless variables this 
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expression reads 

/^o"W 1^2 2 J' 

where h = ho/^/MVoL'^, u"^ = ML'^Uq/Vq. The effective frequency u completely determines 
the classical dynamics. On the other hand, the effective Planck constant h is given by an 
independent combination of parameters. 

One can hardly hope to extract directly the additional factor h}^'^ from the experimental 
data on transmission probability: it is almost impossible to identify the weak power-law 
dependence on top of the leading semiclassical exponent. We suggest an indirect method. 
Namely, consider the quantity 

FQM = -h\og{V/h'/'). 

In the regime of potential tunneling (A oc h^^'^) Fqm is almost independent of K at small 
values of the latter. On the other hand, Fqm — — |^log /?.+ const whenever the new tunneling 
mechanism is involved. The difference between the two cases is seen in Fig. [3l where the 
dependences of Fqm on the total energy E are shown for several values of h. The graphs in 
Fig. [3] coincide at energies somewhat smaller than E^. (say, dX E < 1), while a.t E > E^ a clear 
difference between the graphs appears. We remark that the change in the behavior of the 
exact tunneling probability is gradual, in spite of the fact that the complex trajectories have 
distinct structure aX E < Ec and E > Ec. We discuss this point and derive the appropriate 
uniform formula in Sec. 14.31 

1 

0.8 

^ 0.6 
tin 

0.4 
0.2 


0.95 1 1.05 1.15 1.2 

E 

Figure 3: The probability logarithm Fqm = —h\og{V/h}/'^) plotted as a function of total 
energy for several values of h and Ey = 0.05. Points represent the exact quantum mechanical 
results; the interpolating lines are drawn for convenience. The critical energy is shown by 
dashed vertical line. 
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Another signature of the sphaleron-driven mechanism was first pointed out in Refs. [251 
150] . One notes that the second stage of sphaleron-driven transition, the decay of the 
sphaleron orbit, proceeds classically and does not affect the leading suppression exponent 
F of the probability. In addition, the sphaleron, being unstable, can evolve at the classical 
level into the out-states with different values of oscillator energy E^. Classical trajectories 
corresponding to these evolutions are obtained by adding small momentum in the direction 
of the out-region at different points of the sphaleron orbit. One concludes that in the case of 
sphaleron-driven tunneling the distribution over final oscillator energies is almost constant in 
some region <El < E^2- The latter region corresponds to the decays of the sphaleron 
along different classical trajectories. 

Note that the above feature is in sharp contrast with the properties of final states in the 
standard case of potential tunneling. Namely, in a typical situation the complex trajectory 
describing transmission through the barrier is unique, and the corresponding out-state wave 
function forms sharply peaked Gaussian distribution around some optimal value E^ = (E^). 

To illustrate explicitly the effect of anomalously wide final states in the case of sphaleron- 
driven tunneling, we consider transitions between the exclusive in- and out- states which 
have definite energies of y oscillator, Ey and E^ respectively, and the same total energy E. 
Then, we fix the initial state {E and Ey) and analyze the dependence of the exact exclusive 
probability Ve on E^^. This dependence is shown in Fig. HJin logarithmic scale for several 
values of E. One immediately sees in Fig. that the width of the out-state distribution 
grows as the value of total energy approaches Ec{Ey) from below. In particular, a flat plateau 
gradually develops in the right side of the distribution. At energies higher than critical the 
plateau is wide and corresponds to the maximum probability of tunneling. Moreover, the 
graphs become flatter as the value of h decreases, see Fig. |Dd. 

One sees another feature of the new tunneling mechanism: the short-scale fluctuations 
in the right and left parts of the plateaux in Fig. |Dd. This is the hallmark of quantum 
interference phenomena, which seem to be important for complete understanding of exclusive 
processes at E > Ec{Ey). We discuss this point in Sec. O 

4 Modified semiclassical technique 

In this section we describe the semiclassical technique adapted to the analysis of sphaleron- 
driven tunneling. We start by reviewing the path integral derivation of the standard method 
of complex trajectories [22j. Then, we manipulate with the path integral and obtain the 
modified semiclassical expressions in the case of sphaleron-driven tunneling. 

For simplicity we assume that the system undergoing tunneling transition is similar to the 
model of Sec. [2l Throughout this section we consider tunneling between the asymptotic in- 
and out- regions of two-dimensional waveguide potential, where the in-state of the process 
\E, Ey) is fixed and the final state is inclusive. It is worth noting that both the standard 
and modified semiclassical methods are completely general and the semiclassical formulas 
of this section can be generalized to other systems. In particular, the modified method was 
applied to the case of chaotic tunneling in Ref. [16] and to field theory in Ref. [32] . 
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Figure 4: Distributions of the logarithm of exclusive tunneling probability over the out-state 
quantum number E^^. The graphs are plotted at Ey = 0.05 and: (a) h = 1/40 and four 
values of total energy; (b) E = 1.2 and different values of h. Note that Ec{Ey = 0.05) ~ 1.1. 



4.1 The standard method 

Semiclassical calculations within the method of complex trajectories proceed as follows. One 
reduces the problem of computing the tunneling probability to a problem of finding the com- 
plex trajectory x^'^\t) — complex solution to the classical equations of motion with certain 
boundary conditions. In practice this solution is obtained numerically. Then, tunneling 
probability is given by Eq. ([T]) where F and A are certain functionals of x^^\t). In this 
section we derive the boundary conditions for x^^\t) and expressions for the functionals F, 
A in the standard case of potential tunneling. 

In order to compute the inclusive tunneling probability we first obtain the semiclassical 
expression for the final state of the tunneling process. The wave function of the final 
state has the form, 



^-iH{tf-u)/h 



Xi) ^i{Xi) 



(3) 



where x = {x, y), while \E'i(a;j) = {xi\E, Ey) is the in-state wave function. Below we assume 
implicitly that and ^/ have support in the in- and out- asymptotic regions respectively. 
One uses the path integral representation for the quantum propagator in Eq. ([3]) and writes. 



dxi "^iixi) / [dx 



JS[x]/h 



(4) 



where 5* stands for the classical action of the system. One observes that at small h the 
integrand in Eq. (jlj) contains fast-oscillating exponent; thus, the respective integral can be 
evaluated by the saddle-point method. To keep the discussion short, we defer the details 
of the saddle-point integration to appendix [Aj here we quote the result. One finds the 
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extremum of the leading exponent in Eq. (jl]), which is represented by the trajectory x^^\t) 
going between the in- and out- regions. This trajectory is generically complex. It satisfies 
the classical equations of motion SS/6x{t) = and arrives at a given point x = Xf oXt = tj. 
The boundary conditions at t = for x'^'^\t) are obtained from the saddle-point integration 
over Xi] they fix the values of in-state quantum numbers, 

i?, = (l/,2 + ^V)/2, E = xy2 + Ey, (5) 

where the subscript i marks the quantities evaluated at t = tj. For brevity we omit the 
superscript (s) of the semiclassical trajectory in Eq. ([5]) and in what follows. 

As the result of integration in Eq. (jl]), one finds the semiclassical wave function of the 
final state, 

^f{xf) = D-'/' ■ exp !^^{S[x] + BM) + ^1 , (6) 

where Bi is the in-state contribution to the exponent and D represents the prefactor deter- 
minant, see Eqs. fl50|) . fl58|) in appendix |A] for explicit expressions. Note that the leading 
exponent S + Bi in Eq. is evaluated on the saddle-point trajectory x{t). 

The inclusive probability V of transmission is equal to the total fiuxEl of the out-wave ([6]) 
through the distant line Xf = x^^ , where x^^ is large and positive. Semiclassically, one 
writes, 

|2 



V = J dyf\^f{xf)\'Rexf, (7) 

where we used the fact that dS/dxf = Xf. The integral in the above expression is again 
computed by the saddle-point technique. In appendix |X] we show that the extremum of the 
leading exponent in Eq. ([7]) is achieved when 

Xf = x^f^ , Im?// = Imi/f = . (8) 

Equations ([H]) fix the boundary conditions aXt = tf for the semiclassical trajectory. 

After the saddle-point integration in Eq. ([7j) one finally arrives at the familiar semiclas- 
sical expression ([1]) for the tunneling probability, where the leading exponent is 

Fpot = 2\m{S + Bi) . (9) 

Note that we mark all the standard semiclassical expressions with the subscript pot which 
stands for "potential tunneling" . 

The prefactor Apot is computed as follows (see appendix |X] for the derivation). One 
finds two independent perturbations 5x'^^\t) and 5x'^'^\t) in the background of the complex 
trajectory x{t). These perturbations satisfy the linearized classical equations of motion, 

+ V"{x{t))5x^''^ = , n = 1, 2 (10) 



^We use the in-state with the unit flux normaUzation. 
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with certain Cauchy datc(§ at t = tf. After evolving Sx^'^^t) back in time from t = tj to 
t = ti, one computes the prefactor by the formulc0 

Apot = I = , (11) 

J47rIm(5Ej^[5a;(i)] • 5E*y[5x^^)]) 

where the hnear functional 

5Ey[5x\ = yiSiji + u'^ViSyi (12) 

measures the change in the initial oscillator energy Ey due to the perturbation Sx{t). We 
stress that 5Ey[5x^^'>] involves perturbations in the in-region, while the Cauchy data for 
5x^"'\t) are set at t = tj. We also note that the prefactor ffTTj) is explicitly proportional to 
h}/"^] this fact was used in the previous section. 

The standard semiclassical calculation is summarized as follows. One finds the com- 
plex trajectory x{t) satisfying the classical equations of motion with the boundary condi- 
tions ([5]), ([8]). Our numerical method for finding the trajectory is presented in appendix [Bl 
The suppression exponent Fpot of the probability is given by the value of the functional ([9]) 
on the trajectory x{t). Then, one considers the linear perturbations around the semiclassical 
trajectory and finds the prefactor Apot using the expression (JTTll . 




(a) (b) 

Figure 5: The semiclassical (lines) and exact quantum mechanical (points) results for (a) 
leading suppression exponent and (b) prefactor; Ey = 0.05. Errorbars represent uncertainty 
of the fit (USD. 

Before proceeding to the case of sphaleron-driven tunneling, we demonstrate explicitly 
that the semiclassical expressions ([9]), (|TT1) produce correct values of suppression exponent 

^First, the perturbations are real at t = tf. Second, they do not change the value of total energy, 
i5i?[(5a;("'] = 0. Third, Q{Sx'^^\ Sx''^^) = 1, where Q is the canonical symplectic form. 
^As discussed in appendix |Al this formula is canonically covariant. 
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and prefactor. To this end, we calculate the exact probability of transition by solving nu- 
merically the stationary Schrodinger equation (see Refs. [TB] for the numerical method). 
The exact values of V are computed at several h. Then, the dependence V{h) is fitteq§ 
with the formula 

- h\og{V/W) = Ffu - h\og Afu + h^Cfu , (13) 

where 7 = 1/2 and the last term accounts for the higher-order semiclassical corrections. The 
fit produces the "exact" values Fju, Aju of the suppression exponent and prefactor; they 
should coincide with the corresponding semiclassical quantities. In Fig. [5] we compare the 
semiclassical results computed by Eqs. ffTT]) with those extracted from the fit f[T5]) . One 
observes remarkable agreement. It is worth noting that the fit flTSl) is extremely sensitive 
to the assumed ^-dependence of the prefactor. In particular, if one erroneously uses 7 = 
or 7 = 1 in Eq. (fT3ll . the value of prefactor extracted from the fit becomes close to zero or 
extremely large. Hence, the graph in Fig. [5]d confirms, in particular, the qualitative formula 



4.2 Modification 

At high energies tunneling proceeds by the new mechanism based on qualitatively new prop- 
erties of semiclassical trajectories. Namely, at ii^ > Ec{Ey) the trajectories get attracted to 
the unstable sphaleron orbit and thus become unstable themselves. 

The instability of complex trajectories sets obstacles for the semiclassical description. The 
most important difficulty is related to the calculation of the prefactor Apot- Equation flTTl) 
implies that Ap^ is inversely proportional to the values of linear perturbations 5x^'^'' at t = tj, 
while the Cauchy data for Sx^"^^ are set at t = tf. On the other hand, linear perturbations 
in the background of unstable trajectory contain exponentially growing part. Thus, at 
E > Ec{Ey) when the complex trajectory spends infinite time interval in the vicinity of the 
sphaleron, the formula fllip gives Apot = 0. This means that Eq. fill I) is incorrect in the case 
of sphaleron-driven tunneling and suggests that the prefactor is suppressed by an additional 
power of h. 

The main idea of the modified semiclassical method was proposed in Ref . [28] ; it is close 
in spirit to the constrained instanton technique of Ref. [IQ] . Namely, we evaluate the path 
integral (jlj) for the tunneling amplitude in two steps. First, we integrate over paths spending 
a given time r in the vicinity of the sphaleron. Second, we integrate over r. 

The above manipulations with the path integral lead to the following method. At the 
first step we obtain certain modified boundary value problem for a family of complex tra- 
jectories labeled by the parameter r. These trajectories are stable and interpolate between 
the asymptotic regions x — > ±oo. The second step produces expressions for the suppression 

^To be precise, we use three values h = 1/70, 1/40, 1/30 at E > 0.75, two values h 1/20, 1/10 at 
0.55 < E < 0.75 and only one value ^ = 1/10 at £' < 0.55. This choice is dictated by the limitations of the 
numerical method which does not allow to perform computations when the value of the tunneling probability 
is too small. 

^At E < 0.75 only two values of h were considered, and we set Cfu — 0- At E < 0.55 (one value of h) we 
were unable to extract the prefactor from the quantum mechanical simulation. 
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exponent Fgph and pref actor Asph in the sphaleron-driven case. These expressions relate the 
values of F^ph and A^ph to limits r +00 of certain functionals evaluated on the modified 
trajectories. 

One introduces the functional r = Tj„t[a;] which, roughly speaking, measures the time 
spent by the path x{t) in the region of non-trivial dynamics. We call Tint interaction time. 
It has the following properties. First, Tint is positive-definite for real paths. Second, it is 
finite for any real path satisfying x ±00 as t ^ ±00 and infinite otherwise. The simplest 
choice is 

Ti„tM = j dtf{x{t)), (14) 
where the function f{x) > vanishes at a; ^ ±00. We use 

/(a;)=exp {-{x + yf/2} 

in the model 

Consider the path integral (jlj) for the final state. One inserts into the integrand of Eq. (jlj) 
the unity factor 

f+00 f+oo p-ioo -J 

Jo Jo Jioo ^T^ri 

where the Fourier representation of the 5-function was used in the second equality. Expres- 
sion (jlj) takes the form, 



A{S[x\+ieTir,t[x\)/h 



(16) 



where we changed the order of integrations. One notes that the integral in brackets is exactly 
the same as in Eq. (jlj) up to the substitution 

S[x\ ^ S,[x] = S[x] + ieT,^,[x] . (17) 

This integral is evaluated by the saddle-point method in the same way as the integral in 
Eq. (jl]). Namely, one finds the regularized trajectory x^{t) which extremizes the modified 
action and arrives at the point Xf at t = tf. The initial conditions for the trajectory are 
still given by Eqs. since the evolution in the in- region is not affected by the functional 
Tint- The result of integration in Eq. flTBl) is 

dr j^^ ^e-/^ ■ . |-(SJa;,] + B,[x,]) + -| , (18) 

cf. Eq. on]). The prefactor in this equation is given by the same determinant formula as 
D, but with the substitution S S^, x{t) ^ x^{t). 

Let us remark on the representation (fTSi) . One keeps in mind that the integrand in 
Eq. (fT8|) accounts for the contribution of paths which spend a given time r in the region 
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of finite x (interaction region). In particular, this is true for the saddle-point trajectory 
Xf:{t). The latter interpolates directly between the in- and out- regions x ±00 and thus 
is stable. Note that the stabilization of complex trajectory is achieved by the modification 
of the classical equations of motion. Namely, the substitution f|T7j) modifies the potential of 
the system 

V{x) V{x)-ief{x) . 

We will see below that the relevant values of e are real; thus, x^(t) describes evolution in 
complex potential. 

The rest of the calculation proceeds as follows. One evaluates the saddle-point integral 
with respect to e. The integral over interaction time is kept in front of the formula. This 
ensures stability of complex trajectories. The resulting expression for f is substituted into 
the tunnehng probability A subtle point is that, since V involves the square of the out- 
state, one obtains at this stage two integrals over interaction times r, r', where the latter 
comes from One of these integrals can be computed by the saddle-point technique. 
Indeed, returning to the original expression for the tunneling probability in terms of the 
integral over real paths, one sees that fixing the sum r+ = (r -|- t')/2 is sufficient to make 
both interaction times r and r' finite. Thus, we change the integration variables to r_|_ and 
r_ = r — r' and evaluate the saddle-point integrals over r_ and over the final state. In this 
way we are left with the single integral over r+. 

We leave the details of the above computation to appendix [C] and discuss the result. 
First, one arrives at the saddle-point conditions 

ReT,nt[x,] = T+ , e = e*, (19) 

which come from the integrals over e and r_ respectively. The integral over out-states 
produces, as before, the boundary conditions ([8]) at t = tf for x^{t). Note that the first of 
Eqs. din]) implies, in particular, that x^{t) is stable. The result for the probability is 



r 



'TT 



h 



de 



1/2 

■ ApoM e-(''-*'^-'^"+^/\ (20) 



where the suppression exponent Fpot,e and prefactor Apot^e are computed by the same for- 
mulas ([9]) and (fTTj) as before, but with the substitution S — > S^. Note that the latter 
substitution implies that both the classical equations of motion and linearized equations 
(ITU]) are modified. 

We now proceed to the second step of the calculation and consider the integral over the 
interaction time r_|_. One makes an important observation: the values of r_|_ and e(r_|_) 
are related by the Legendre transformation. Indeed, by construction the configuration 
{Xe{t), e(T+)} corresponds to the extremum of the leading exponent = Fpot,e — 2er+ 
in Eq. (|20|) . and the respective derivatives of F^ are equal to zero. Thus, 

£ = ^(^ioM - 2.r,) = -2. . (21) 
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where only the exphcit dependence of on r_|_ was taken into account in the last equal- 
ity. Due to the property (ETj) . the integral in Eq. (1201) is saturated at e = 0. This point 
corresponds to the original semiclassical equations: recall that the modification term in 
the classical action, Eq. f|T71) . is proportional to e. One concludes that the integral for the 
tunneling probability is saturated in the vicinity of the original complex trajectory at e = 0. 

So far in our calculation we did not make any reference to the particular tunneling 
mechanism. Thus, Eq. fl2U]) can be used in cases of both potential and sphaleron-driven 
tunneling. The difference between the two mechanisms becomes crucial in the evaluation of 
the integral over r_|_. In the standard case of potential tunneling the trajectory at e = is 
stable and corresponds to the finite value of interaction time r+; one takes the integral in 
Eq. (120|) by the saddle-point method and arrives at the expressions ([9]), (fTTj) from the previous 
subsection. The case of sphaleron-driven tunneling is considerably different, because the 
time interval spent by the trajectory in the vicinity of the sphaleron tends to infinity as 
e —>■ +0. Thus, the integral in Eq. fl20l) is saturated by the end-point of the integration 
interval r+ — > +oo. Using the appropriate asymptotic expressioio for the integral, one 
obtains Eq. ([1]) with 

Fsph = lim , (22a) 
Asph = n'/^ hm ^ , (22b) 



dt 



where we mark the quantities corresponding to the new mechanism with the subscript sph. 
Note that the prefactor Apot,e is computed by the formula ffTTl) with modification ([T7|l . while 
the exponent 

Fe = Fpot,e - 2er+ = 2lm{S[x,] + B,[x,]) 

is given by the value of the original action on the modified trajectory. Let us remark that 
the limit e — > +0 in Eqs. (1221) does actually exist; this is shown analytically in appendix [Dl 
One observes that the expression (122bp for the prefactor is very different from that in the 
case of potential tunneling. In particular, Agph oc Jl^^'^Apot- 

To summarize, we derived the following method of calculating the probability of sphaleron- 
driven tunneling. One modifies the classical action of the system by adding purely imaginary 
term proportional to the small regularization parameter e > 0, see Eq. f[T7l) . Then one solves 
the modified equations of motion with the original boundary conditions and finds 

the modified complex trajectory x^{t). This trajectory interpolates between the asymptotic 
regions x ±oo and is stable. The modified values of the suppression exponent Fpot,e and 
prefactor Apot,e are computed by the same formulas as before, Eqs. ([9]) and (ITT]) , but with the 



^°This expression is derived as follows. One moves the leading exponent in Eq. ([20]) under the differential 
using the relation 2e • exp{— Fe/fi,} (iT+ = hdcxp{—Fg/h} and integrates by parts. After integration the 
leading semiclassical approximation is given by the boundary term at t+ —f +oo] the boundary term at 
r+ = and the remaining integral over r+ are exponentially and power-law suppressed respectively. 
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substitution S ^ S^. The final result for the tunneling probability is obtainec^ in the limit 
e — s> +0, see Eqs. (1221) . We call the above modified semiclassical method by e-regularization 
technique 
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Figure 6: (a) Suppression exponent and (b) prefactor in the case of sphaleron-driven tun- 
neling, E > Ec{Ey)] Ey = 0.05. The points are extracted from the fit f|T3|) . while the lines 
stand for the semiclassical results, Eqs. (l22l) . The vertical dotted lines bound the range of 
energies for sphaleron-driven tunneling. 

We perform straightforward check of the modified semiclassical method by comparing 
the semiclassical predictions (122]) with the results of the exact quantum mechanical compu- 
tations. The latter are used to extract the values of the suppression exponent and prefactor 
by the fitting procedure described in the previous subsection, where 7 = 1 in Eq. f|T3|) . 
The comparison is shown in Fig. [61 The observed agreement between the semiclassical and 
quantum mechanical results supports the modified semiclassical technique. In particular, 
we checked that the fit ffTSl) produces unacceptably large values of the prefactor if one er- 
roneously assumes the same fi-dependence 7 = 1/2 as in the case of potential tunneling. 
Thus, the scaling Asph oc ^ is confirmed. 



4.3 Uniform approximation 

Our expressions for Apot and A^ph imply apparent discontinuity of the semiclassical probabil- 
ity across the critical energy; on the other hand, the exact quantum probability is a smooth 
function of energy. As a consequence, the fl-dependences Apot oc fi^^^ and A^ph oc h fail to 
describe the quantum mechanical data in the immediate vicinity of Ec{Ey). [This is seen in 
Figs. [5]d, [6]d, where the quality of the fit (fT3!) becomes worse as E ^ Ec{Ey).] One observes 
that both the standard and modified formulas are invalid at ~ Ec{Ey). 

^^In practice the limit in Eqs. ([22]) is taken by considering small values of the regularization parameter, 
£ ^ 10^^. At these e the values of the suppression exponent Fgph and prefactor Agph stabilize at the level of 
accuracy 10~^. 
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In this section we derive the uniform asymptotic formula for the tunnehng probabihty 
which is apphcable in the vicinity of Ec{Ey). Our formula has the form (cf. Ref. [H]), 



um 



Mpot ■ Vpot at E < EciEy) , 
Msph ■ 'Psph at E > Ec{Ey) , 



(23) 



where Aipot and Aisph are the correction factors in the cases of potential and sphaleron- 
driven tunneling respectively. We will find that A4pot,sph ~ 1 at — Ed ^ Jl^^'^] thus, the 
formula fl23l) is relevant in the small region of width AE ~ h}^'^ around the critical point. 
We stress that the uniform probability Vuni is continuous at E = Ec{Ey). 

We obtain the desired approximation by examining the integral over r+ for the tunneling 
probability, Eq. fl2Ul) . Recall that Eq. fl2U]) is applicable in both cases of potential and 
sphaleron-driven tunneling. To make the discussion transparent, we change the integration 
variable to 



Note that the limiting values r+ = and r_|_ +oo correspond respectively to w = Wq > 
and w — > +0. In new terms the integral fl20p takes a particularly simple form. 



where the leading semiclassical exponent is now considered as function of w. 

The difference between the two mechanisms of tunneling is now understood as follows. 
At small energies the integral fl25l) is saturated by the saddle point w = Wg > 0. The value 
of Wg decreases with energy, so that at E = Ec{Ey) the saddle point crosses the boundary 
w = and leaves the integration interval. At ii^ > Ec{Ey) the saddle point Ws is situated 
outside the region of integration. 

Consider the Taylor series expansions of the semiclassical exponent around the points 
w = Ws and w = 0, 



where the primes denote derivatives with respect to w and we marked by Fpot, Fgph the values 
of the exponent at w = Ws and w = 0. The semiclassical expressions of Sees. 14.11 and 14.21 
are obtained from the expansions 0261) and (127|) respectively. Namely, at E < Ec{Ey) one 
implements the saddle-point method, i.e. substitutes Eq. fl26l) into Eq. fl25l) and extends the 
interval of integration to the entire w axis. At energies higher than critical the minimum 
value of the exponent F^ is achieved at w = 0, and one uses the expansion fl27j) . where only 
the zeroth- and first-order terms are kept. It is straightforward to check that in this way one 
obtains the expressions ©, (HH) of Sec. O and of Sec. S^i 

One observes that the above two integration methods are not applicable if the saddle point 
Ws is close to the end-point w = 0. Indeed, in the saddle-point technique at E < Ec{Ey) the 




(24) 




(25) 



F,{w) = Fpot + F"{ws) ■ (w - Wsf/2 + 0{{w - Wsf) , 
FM = Fsph + F'iQ) ■ w + F"(0) • w^/2 + 0(«;3) , 



(26) 
(27) 
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interval of integration cannot be extended to the entire w axis, since the contribution from 
the additional interval w < is not negligible. In the end-point integration at E > Ec{Ey) 
the third term in Eq. fl27p is not small in comparison with the second term, because the 
first derivative F'{0) vanishes in the limit — > 0. Given these reasons, one easily remedies 
the formulas keeping the finite integration interval at E < Ec{Ey) and three terms in the 
expansion fl271) at energies higher than critical. The resulting expressions for the correction 
factors are 

Mpot = ^{l + ^{>ipot)} , where >Cpot = Ws\/ F"{ws) /2h , (28a) 

Msph = V^y^sph {1 - H>^sph)} ■ , where Xsph = F'{0)/^/2hF"{0) . (28b) 

Here $ is the Fresnel integral, ^{z) = 2/y^ dte~^^. 

It is straightforward to check that the factors ( l28l) have the required properties. First, the 
uniform formula (123|) is continuous at = Ec{Ey) by construction. Indeed, in this case Ws 
and F'{0) are equal to zero and the expressions fl26l) . (!27|) used in the integration coincide. 
Second, at \E — Ec{Ey)\ ^ h}^"^ the arguments of the Fresnel integrals are large. Using 
the asymptotics of ^{z), one finds that M.pot,sph ~ 1 outside the immediate vicinity of the 
critical energy. Third, one notes that in the region \E - Ec{Ey)\ ~ n}/^ the "potential" and 
"sphaleron-driven" parts of the uniform formula coincide up to higher-order semiclassical 
corrections. Indeed, in this region the central points w = Wg and w = of respective Taylor 
expansions are parametrically close to each other, Wg ~ h^^"^] thus, the results obtained from 
Eqs. ( l26|) and ( !27|) are close as well. 

Numerically, one extracts the quantities entering the correction factors ( l28l) by studying 
the dependence of F^ and Apot,e on e. We discuss this calculation in appendix [El 

To summarize, we derived the continuous asymptotic formula for the tunneling proba- 
bility, Eq. (l23l) . which works at energies close to critical and interpolates between the two 
semiclassical expressions corresponding to the cases of potential and sphaleron-driven tunnel- 
ing. In Fig. [7] we compare the uniform approximation (l23l) (solid line) with the semiclassical 
probabilities Vpot and Vsph (dashed lines), as well as with the exact quantum probability 
(points). 

5 Exclusive processes 

Here we study semiclassically the effect of the new tunneling mechanism on exclusive pro- 
cesses, i.e. processes with completely fixed out-states. We discuss the apphcation of the 
modified semiclassical technique to the exclusive case and obtain expressions, analogous to 
Eqs. f l22|) . for the suppression exponent and prefactor of exclusive probability. We show that 
in the semiclassical limit of vanishingly small h the exclusive prefactor A^^sph is proportional 
to h"^ in the sphaleron-driven case. This should be compared with the dependence A^^pot oc h 
in the case of potential tunneling. 
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Figure 7: Uniform approximation for the tunneling probability (solid line) shown at ^ = 1/30, 
Ey = 0.05 in the vicinity of the critical energy Ec{Ey) ^ 1.1. We also plot the semiclassical 
results of Sees. 14.11 and 14.21 (dashed lines) and exact quantum probability (points). The 
critical energy is marked with the vertical dotted line. 

5.1 Exclusive trajectories 

We consider tunneling transitions between the exclusive states \E, Ey) and \E, E^^) specified 
by the same value of total energy E and definite energies Ey, E^ of y-oscillator. The standard 
semiclassical method in the case of exclusive transitions is formulated in Ref. [22]. Its 
derivation is completely analogous to that carried out in Sec. 14.11 for inclusive processes. 
Fixation of the out-state changes the final boundary conditions for the complex trajectory: 
instead of Eqs. ([HD one has, 



The initial conditions remain the same, Eqs. (jSj). The exclusive suppression exponent is 
given by the action functional 



computed on the trajectory, cf. Eq. iQ. Note that the new term Bf in Eq. fl30|) is related to 
the out-state of the process; it is given by the same expression as Bi, but ai t = tf and with 
the out-state quantum numbers E, E^. We do not write here the formula for the prefactor 
Ae,pou it can be found in Ref. [22] • Importantly, this formula implies that A^^pot oc h. 

We apply the above method in the case of potential tunneling, E < Ec{Ey). In Fig. [8^ 
(lines) we plot the out-state distributions of the exclusive exponent (130|) for several values 



Xf = X 



,(0) 



Ei; = iyj + u;Vf)/2. 



(29) 



= 2lm{S[x] + B,[x] - Bf[x]) 



(30) 
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Figure 8: Exclusive suppression exponent Ff. in the cases of (a) potential and (b) sphaleron- 
driven tunneling; Ey = 0.05. The semiclassical results for F^. (lines) are compared with the 
exact data (points). Errorbars represent the uncertainty of the fit ( fT3l) . 



of energy E. The exact results (points) are extracted from the fit (|T3l) with 7 = 1. The 
semiclassical and exact data coincide. 

One observes that well below the critical energy Ec{Ey) ^ 1.1 the function Fe^pot{Ep has 
a clear minimum corresponding to a sharp maximum of the quantum probability. As the 
energy tends to Ec{Ey), a fiat plateau develops in the right side of the graph. As discussed 
in Sec. El this behavior is copied by the exact quantum probability, cf. Fig. H^. 

Before introducing the semiclassical method for exclusive tunneling in the sphaleron- 
driven case, we preview the result for the suppression exponent F^^sph in Fig. [8)d (solid line). 
At E > Ec{Ey) the exclusive exponent is exactly constant in the region Ey^^ < Ef^ < Ey 2] 
clearly, this feature corresponds to a wide and fiat maximum of quantum probability, cf. the 
exact graphs in Fig. |Dd. Thus, the distribution of the exclusive probability over the out- 
state quantum numbers becomes anomalously wide when the sphaleron-driven mechanism 
is involved. So far the semiclassical study of this property was restricted to one-dimensional 
systems with non-autonomous potentials [251 130] . Here we find the same effect in the two- 
dimensional setup of Sec. [2l 

Returning to the semiclassical description of exclusive tunneling processes, we find the 
following manifestation of the sphaleron-driven mechanism. In contrast to the case of po- 
tential tunneling where the exclusive trajectory is unique, dX E > Ec{Ey) there is an infinite 
sequence of complex trajectories corresponding to the same final oscillator energy E'^. In 
Fig. [9lwe plot the first four trajectories for Ej^ = 0.6. One observes the following behavior: 
the trajectories reach the unstable periodic orbit (sphaleron), perform several oscillations 
there (i.e. around the point a; = 0) and then slide off describing the sphaleron decay into 
the final state. Importantly, in order to arrive into the out-state with given E^, the trajec- 
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Figure 9: Four exclusive trajectories at = 0.6. The initial parts of the trajectories are 
indistinguishable on the graph. The in-state quantum numbers are E = 1.2, Ey = 0.05. 

tory must leavj^ the sphaleron at a particular oscillation phase cp. More precisely, there 
are two choice^ of phase per sphaleron period. We conclude that the interaction time r_|_ 
spent by the exclusive tunneling trajectories at given E^ is restricted to two values plus an 
integer number of sphaleron periods. This gives rise to an infinite family of tunneling trajec- 
tories which describe the same exclusive process but differ by the number of "half-period" 
oscillations on top of the unstable periodic orbit. 

To investigate the properties of exclusive trajectories in the sphaleron-driven case, we 
proceed as follows. For each trajectory we compute the value of the interaction time func- 
tional r+ = lie Tint[x], Eq. f|T^ . In this way we characterize the trajectories by points in the 
plane (i?^,r+), see Fig. [TUl 

A comment is in order. The exclusive tunneling trajectories are stable even in the 
sphaleron-driven case. Thus, they can be found without e-regularization. Still, as we will 
discuss shortly, it is convenient to use the modified semiclassical technique at the intermedi- 
ate steps of the computation and remove the regularization afterwards. To avoid confusion, 
let us stress that the solid line in Fig. [10] corresponds to trajectories which are obtained 
a/i(er removal of the regularization. Consequently, the functional Tj„t[a;] does not enter the 
equations of motion for these trajectories and is used only to characterize their temporal 
behavior. 

From Fig.[TO]one sees that the exclusive trajectories are naturally divided into two classes. 
The trajectories from the first class lie in the interval Ey^ < E^ < Ey2 corresponding to the 
plateau in Fig. In Fig. (TU] they form a sin-like curve extended to the infinite values of r+. 
All these trajectories describe creation and subsequent decay of the sphaleron. Moreover, 
we find that the latter decay proceeds classically since the imaginary part of the trajectories 

^^This notion can be given precise meaning by saying that the trajectory leaves the sphaleron once the 
distance \x{t) — Xsph{t)\ between the trajectory and the sphaleron orbit reaches a certain value S <S^1. 
^•^This follows from the fact that the final oscillator energy is a periodic function of ip; thus, equation 
if ) = const has (at least) two solutions. 
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Figure 10: The curve in {E^^, r+)-plane (sohd hne) representing exclusive tunneling trajec- 
tories. Circles correspond to the trajectories from Fig. O The dashed line represents the 
modified inclusive trajectories. It is almost coincident with the sin-like part of the solid line. 
Arrows indicate the direction of decreasing e. The in-state quantum numbers are the same 
as in Fig. [91 

becomes small after one sphaleron oscillation. As a consequence, the value of the functional 
( I30ll is almost independent of the individual trajectory from the first class. Besides, it is 
clear from the figure that these trajectories form an infinite sequence of branches marked 
with the integer number m of "half-period" oscillations in the vicinity of the sphaleron. 

The trajectories from the second class represent the "wings" E^ < Ey^^, E^^ > Ey ^ of 
the out-state distribution in Fig. [Hb- They correspond to the case when the decay of the 
sphaleron orbit into the out-state with given Ej^ cannot proceed classically. Consequently, 
the probabihty of this decay is exponentially suppressed. Due to the additional suppression, 
the exponent (l30l) strongly depends on the out-state at Ej^ < Ey^^, E^ > Ey2- Note, 
however, that the sphaleron still serves as the mediator of the two-stage tunneling process; 
hence, the second-class trajectories with fixed E^ form an infinite sequence marked with the 

topological number m, see Fig. [TDl The values of the suppression exponents -Fijof calculated 
on trajectories with different topology and given E^ are almost degenerate. 

In practice the exclusive trajectories are conveniently found using the regularization 
method of Sec. 14.21 The procedure is based on the following observation. Consider e- 
regularized trajectories corresponding to the inclusive tunneling process. They describe 
creation and subsequent classical decay of the sphaleron. The final state of the decay de- 
pends on the value of e. Changing e one covers the whole range of final oscillator energies 
Ey^^ < E^ < Ey 2 accessible in the classical sphaleron decay. This consideration is illustrated 
by the dashed curve in Fig. [10] which represents the modified inclusive trajectories at dif- 
ferent values of e in the {E^, r_|_) - plane. We see that the graph closely follows the sin-like 
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curve of exclusive trajectories from the first class, and the value of e decreases towards large 
r_|_ (along arrows). Thus, the modified solutions with different e form a single branch which 
smoothly interpolates between the branches of exclusive trajectories. 

Numerically, we exploit the above property by applying the deformation procedure of 
appendix [HI Namely, we start with the modified trajectory at a given Ej^ = Ey ^. Suppose it 
has topology m. Then, the trajectory with topology m + 1 (m — 1) is obtained by decreasing 
(increasing) the value of e until the final oscillator energy arrives to Ey^ again (see Fig. [TU]) . 

Repeating this procedure, we find the sequence of modified trajectories at E^ = Ey ^. Finally, 
we impose the boundary conditions fl29|) and set e = 0. In this way we find all exclusive 
trajectories from the first class sorted by the topological number m. The solutions at the 
"wings" are obtained by taking the trajectories corresponding to E^^ = Ey^^ i^y,2) ^^"^ 
deforming them by decreasing (increasing) Ej^. 

5.2 Exclusive probability 

Let us derive the expression of the form ([T]) for the exclusive tunneling probability in the 
sphaleron-driven case. We start with the semiclassical formula 

oo 

p. = E-4ii (31) 

m=0 

where the sum runs over all complex trajectories describing the same process. Note that 
the terms due to interference between different trajectories are neglected in Eq. fl3Tl) : we will 
discuss them later. One recalls that the number m of exclusive trajectory increases with the 
time interval r+ spent by the trajectory in the vicinity of the sphaleron orbit. In accordance 
with the new tunneling mechanism this implies that the sum in Eq. flHT]) is saturated at 
m —>■ +oo: the individual suppressions F^^l^ decrease with m and reach the minimum at 
m —>■ +00. This minimum is the overall suppression exponent of the process, 

Fe,,p,= hm (32) 

Note that the value of F(,^sph is the same for all trajectories from the first clas^ and equal 
to the suppression Fgph of inclusive tunneling probability. This property gives rise to the 
plateau in the dependence Fe^sphiE^^) in Fig. [8)3. 

Let us now turn to the prefactor. Since the suppressions change at large m in small 
steps, one may be tempted to replace the sum in Eq. (1311) by the integral and evaluate it 
in a straightforward way. However, this replacement is in general incorrect: even for small 
change of ivpoi the change in the exponent e~ ^'P"*' can be large. 

^^One proves this by noting that exclusive trajectories at large m are close to the respective modified 
trajectories, and the limit m —t +00 in Eq. (|32p can be substituted with r+ —t +00. Since the modified 
trajectories sweep the interval Ey ^ < < Ey 2 as T-f grows, the limiting value Fe^sph does not depend on 
EI within this interval. 
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We proceed carefully. In what follows we restrict our attention to the plateau case 
E^^^^ < < Ey^2- starts by relating the limit m +00 of exclusive quantities to Asph, 



lim 4^™) 
nrn /±e,pot 

m— >+oo 



"-^ e,pot 



dEl 



ujA 



sph 



(33) 



This formula is obtained as follows. One changes the integration variables from r_|_ to E^ in 
the expression fl20|) for inclusive probability, 



dEl 



where the derivative is taken along the dashed line in Fig. [TOl Comparing the resulting 
integral with the relation 

dEl 



V 



between the inclusive and exclusive probabilities, one expresses the modified suppression 
exponent and pref actor in terms of -Fg.pot) ^i'^L; 



^(m) 
e,pot 5 
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pot,e 



e,pot 



dr. 



^/h[ de _ 
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dEf 
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(34) 



Note that the value of e in these formulas is fixed by the specification of the final oscillator 
energy E^ and topological number m of the respective trajectory. Now, one notes that the 
limit in Eq. fl22bp can be computed by considering the subclass of modified trajectories with 
fixed E^j. These are close to the respective exclusive solutions; one uses the latter in the 
r.h.s. of Eq. (]22b|) and substitutes the limit e +0 with m — * +00. Then, Eqs. and 
the Legendre transformation (12T|) imply Eq. (1331) . 

Now, we exploit the dependence of the individual suppressions F^'^l^ on m at large m. It 
is shown in appendix [D] that the suppressions approach the limiting value Fgph exponentially. 



^ e,pot sph 



aodd{El) e 



m = 2n 
m = 2n + 1 



(35) 



where the coefficient (3 = PT^ph is related to the positive Lyapunov exponent (3 and pe- 
riod Tgph of the sphaleron orbit. Clearly, P does not depend on the final oscillator energy. 
Substituting Eqs. (l35l) . (133|1 into the formula (!3T|) . we find. 
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(36) 
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Let us concentrate on the first term in braces, the second term is treated in the same way. 
The sum is saturated near the point corresponding to the maximum of the exponent, 



1 h 

no = --5 log' 



(37) 



Generically, no is not integer. Factoring out the value of the summand at n = no, one writes, 
00 »^ 00 

n=[} 



J2 exp (-/3(n - no) - e-^^"""")) 



(38) 



where in the r.h.s. we extended the sum to all integer n by noting that the terms at n < 
are negligibly small. Substituting this relation into Eq. ( 136|) . one finally obtains expression 
for the exclusive prefactor. 
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where 



s{P,0^ J2 exp(-/3n-e 



) 



(39) 



(40) 



is a periodic function of ^ with period p. 

Let us discuss our result. The dependence of the exclusive prefactor (l39l) on h is different 
in the cases /3 <C 1, /3 ~ 1, /? 3> 1. At /5 -C 1 the sum in Eq. (l40l) can be replaced by 
the integral and one obtains s(/?,^) ~ 1//?. Then the ^-dependence of Ae^sph reduces to the 
simple proportionality la"ft0 A^^sph oc h?. In the generic case /? ~ 1 one observes, besides the 
overall scaling A^^sph oc Ji^, the modulation of the prefactor by the periodic function of log^. 
The latter modulation is elusive, however, in models with /? 3> 1. Namely, the periodic 
nature of s{P,^) becomes apparent only at |,^| = | \og{h/aeven)\ ~ P which corresponds to 
exponentially small values of h ^ e~^. 

Realistically, at large /? one works in the regime /? ^ |^|. In this case the formulas ( !32l) . 
( I39l) are not applicable, since they are derived under the assumption no > 1, see Eq. ( 1371) . At 
no ^ 1 one uses the original expression ( l3Til . where the sums over even/odd m are saturated 
by the first terms. Then A^^sph oc h. Let us roughly estimate the relative size of the terms 
with m = and m = 1. One takes 



^ e,pot sph 



^ ^{-^e,pot Fsph) 



where c is a coefficient of order 1, and uses Eq. ( I33l) to estimate the prefactors. This yields 
that the trajectory with m = 1 is relevant when 



^ ^ ^e,pot ^sph 

cP 



(41) 



1^ Recall that A 
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Figure 11: Exclusive suppression exponents in the sphaleron-driven case: the first exponent 
F'e^pot (liiis), the hmiting value Fe^sph (dashed line) and the fit of the exact data with H > 1/120 
(errorbars). The in-state quantum numbers are the same as in Figs. [8}d,|4)d. 



and is negligible at larger h. One concludes that, depending on the value of H, the term with 
m = or m = 1 dominates. 

The characteristic values of the parameter (3 are related to the properties of the unstable 
periodic orbits, which are fixed in the model under consideration. As estimated in appendix 
ID| in the setup ([2]) /5 ~ 24. On the other hand, the numerical quantum mechanical compu- 
tations are feasible only down to h> 10~^. Thus, we are in the regime h S> e~^, where the 
exclusive probability is saturated by the trajectories with m = 0, 1. 

To see this explicitly, we compare the first suppression exponent -Fipo^ and the limiting 
exponent^ F^^sph with the exact suppression. Since in our case A^^sph oc h, one extracts the 
exact suppression from the fit (fT3|) with 7 = 1. We consider separately the exact data with 
h > 1/120 and h < 1/120. Using the first set of data, we obtain the points in Fig. [TTl The 
results of the fit closely follow -Fg and differ substantially from the limiting exponent F^^sph 
in the left part of the graph. One concludes that at h > 1/120 the trajectory with m = 
saturates the tunneling probability. 

Second, we analyze the exact quantum data with h < 1/120. Consider Fig. |3)d where 
the logarithm of the exact tunneling probability is plotted at several values of h. One sees 
that the graph at h = 1/140 is notably different from the graphs at larger h. We attribute 
this difference to the contribution of the trajectory with m = 1. Indeed, at ^ E^^ the 

difference in the suppressions is large, -Ff pot — -Fg.sph ~ 0.2. Then, the estimate fllTl) implies 
that at ^ < 1/120 the first odd trajectory enters into the game. The results of the fit of 
exact data with small h are shown in Fig. [8]d. As expected, they coincide with Fe^sph in the 
leftmost part of the graph. 

Finally, let us briefly discuss interference between exclusive trajectories. Consider first 
the case /? ~ 1. The sum fl40l) for the prefactor is then saturated by the fixed number of 

^^In our case Fe^pot = Fe,sph with good accuracy. 
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terms, An = 0{fP). Since each term corresponds to the complex trajectory, the number of 
trajectories giving substantial contribution into the probability is finite, and the interference 
between the trajectories is important. This gives rise to oscillations in the dependence of 
the probability on the in- and out-state quantum numbers Ey, E^. The period of these 
oscillations tends to zero a.s h —>■ 0; thus, they become indiscernible in the semi classical limit. 
However, at finite h the interference is important. 

In our case of large (3 the exclusive probability is dominated by two complex trajectories, 
and the interference picture is seen whenever the contributions of these trajectories are 
comparable. In accordance with the above discussion, this happens at ^ ~ 1/140 and 
E^ -E^i, see Fig. HJd. The small-scale oscillations in the right part of the plateau in 
Fig. IDd are explained as follows. Let us take a look at Fig. [TOl One observes that at 
-E^ ~ Ey2 the trajectory with m = 1 is almost coincident with the dominant one. Thus, this 
trajectory gives substantial contribution into the probability in the vicinity of Ey2, and the 
interference between the two trajectories is seen in this region. 

6 The limit of small quantum numbers 

According to the common lore low-lying quantum states are "not semiclassical." Indeed, one 
cannot use the semiclassical expressions for the wave functions of these states in the majority 
of applications: at ~ /i the momentum is parametrically small and the semiclassical 
approximation is not justified. One finds, however, that tunneling processes are very special 
in regard of low-lying states. Namely, the semiclassical tunneling probability depends only 
on exponentially small tails of in- and out-state wave functions; these tails can be computed 
semiclassically even at small values of respective quantum numbers. 

In this section we generalize the semiclassical method to the case of tunneling from the 
low-lying in-states of ?/-oscillator, Ey ~ h. At the same time the total energy is assumed to 
be semiclassically large, ~ 1. To be concrete, we take Ey = hw/2, which corresponds to 
the oscillator ground state. Note, however, that the method of this section can be used for 
other low-lying oscillator states as well. 

Let us address the following questions: 

(i) Is it legitimate to use the semiclassical approximation for the wave function of the oscil- 
lator deep inside the classically forbidden region, \y\ ^ h^^"^ ? 

(ii) Is the integral over initial states in Eq. (jl]) saturated deep inside the classically forbidden 
region at small Eyl 

If the answers to the above questions are positive, one can use the semiclassical expressions 
for the probability of tunneling from the ground state of y-oscillator (e.g. Eqs. ([2]), ffTTl) in 
the case of potential tunneling). 

To answer the first question, we compare the semiclassical and exact oscillator wave 
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functions in the case of ground state, Ey = huj/2: 



^ 1/2 



V'».o(a) = ■ exp (-^!/V2ft) . (43) 

In the above expressions \y\ is large and Py,i{y) = ^/2Ey — uP'y'^. Equations (H2l) and (H3l) 
look quite different: the exact wave function involves the factor hr^l^ which is not present 
in the semiclassical expression. However, substituting Ey = hw/2 into the leading exponent 
of Eq. (USD, one finds, 

/y 
dy'Piiy') = Wl2 - ih/A - ih/2 log (2y^/^^ + 0{n^) . 

Thus, up to high-order semiclassical corrections 



TT 



1/4 



'^vM = - --^vM „ _ • (44) 



One concludes that the two wave functions are related by the simple renormalization factor 
(7r/e)V4. 

The relation (jUj) is not surprising. Indeed, the standard derivation of the semiclassical 
wave function (H2|) proceeds in two steps. First, one solves the Schrodinger equation with 
'4'y,s{y) = C ■ e*°"'^^)/'^ considering a/h^ 1. This is certainly valid deep inside the classically 
forbidden region, even for Ey = hjj/2. The second step is the evaluation of the normalization 
constant C by taking the integral / \tpy^s{y')\'^dy' = 1. At small Ey the latter integral is 
saturated at y ~ 0, i.e. right in the vicinity of the turning points, where the semiclassical 
expression P2|) is not applicable. Consequently, the semiclassical calculation produces an 
incorrect value for the constant C a,t Ey = fiw /2. Equation (jUj) shows that the correct value 
is (tt/c)^/^ times larger than the one obtained semiclassically. 

We have the following answer to the question (i): the semiclassical expression flTU]) can 
be used at Ey = huj/2 deep inside the classically forbidden region; however, the final result 
for the probability should be multiplied by the correction factor (vr/e)^/^. 

Let us address the question (ii). Consider the complex trajectory x{t) in the in-region. 
One finds, 

x{t) ^ p.^^i{t - U) + Xi , ?/(t) ^ ae"*"^* + ae*"^* as t ^ -00 . (45) 

The initial boundary conditions ([5]) guarantee that the quantities px,i = \/2{E — Ey) and 
aa = Ey/2uj'^ are real. Therefore, one can define two real parameters T, 9 by the relations 

Imxi = -Px,iT , a* = ae-^'^^'^ . (46) 
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As discussed in Refs. [381 126], these parameters are in one-to-one correspondence with the 
in-state quantum numbers E, Ey. In other words, T and 6 provide an alternative parameteri- 
zation of tunnehng trajectories. Note that T = 6 = represent classically allowed transitions, 
x{t) G M. On the other hand, the limit 6 — > +oo corresponds to Ey —>■ 0. Indeed, in this 
limit one obtains a —>■ and a finite [38l [26] , which are the Feynman boundary conditions 
for tunneling from the ground state. From Eq. fj^5|) one finds that \yi\ \d\ ^ h}^"^. Thus, 
the integral over initial states in Eq. (jl]) is saturated deep inside the classically forbidden 
region, where the semiclassical expression for the in-state wave function is trustworthy. 

One concludes that, apart from the additional multiplier (vr/e)^/^, the semiclassical ex- 
pressions for the tunnehng probabihty ([1]), such as Eqs. Qj, (fTTl) . are still applicable at 
Ey = nw/2. 

Note that in the considered case of tunneling from the ground state the expressions ([9]), 
(ITT]) depend on h in non-trivial way through Ey = huj/2. It is convenient to extract this 
dependence explicitly and bring the expression for the tunneling probability into the form 
([1]) with F independent of h and A having only the power-law dependence. This is done in 
appendix [Fl the result is 

Fpotfi= lim Fpot, Apotfl = {'K/Kf^ e'^'^ lim A^ot ■ (47) 

where Fpot and Apot are the standard semiclassical expressions for the suppression exponent 
and prefactor. The quantity entering Eq. W7\i is extracted from the smal\-Ey asymptotic 
of the leading exponent 

Fpot = Fpotfl + ^ log(2^,/^) - ^^Ey + 0{El) . (48) 

U] UJ 

Let us remark on Eqs. (147]) . First, note that Ap^tfi contains the additional factor h^^^"^ 
as compared to the case of highly excited in-states. Second, we did not use the dynamical 
properties of complex trajectories in the derivation of Eqs. ( HTj) . Thus, the above expressions 
are valid both for inclusive and exclusive processes. They also hold in the case of sphaleron- 
driven tunneling, where one substitutes Fpot Fgph, Apot Asph in Eqs. fHTj) . In particular, 
for the prefactors of inclusive processes one has Apotfl oc hP and Asph,o oc h^^'^. 

Finally, it is worth mentioning that the first of Eqs. fHTj) . namely, the limiting relation 
between the suppression exponents of tunneling from the low-lying and highly excited in- 
states is known in field theory as the Rubakov-Son-Tinyakov conjecture [M]- We proved 
this conjecture in quantum mechanical setup. 

We close this section by comparing the semiclassical results for the suppression exponent 
and prefactor, Eqs. (HTj) . with the results extracted by the fit Eq. (fT3|l from the solution of 
the Schrodinger equation. The comparison in the cases of potential and sphaleron-driven 
tunneling is presented in Figs. ITIk . ITIb for inclusive and in Figs. [T^ . [T5b for exclusive 
processes. In the latter case we compare the exact suppression exponent with the suppression 
of the first exclusive trajectory. One observes nice agreement. 
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Figure 12: Comparison between the exact (points) and semiclassical (lines) results for the 
suppression exponent and prefactor at Ey = huj/2 in the cases of (a) potential and (b) 
sphaleron-driven tunneling. Note that Ec{0) ~ 1.1. Errorbars represent uncertainty of the 

fit (USD. 





Figure 13: Comparison between the exact (points) and semiclassical (lines) results for the 
suppression exponent of exclusive processes at Ey = fuj /2 in the cases of (a) potential and (b) 
sphaleron-driven tunneling. In the case (b) the exact data coincide with the first exponent 
Fgpgj Q. Errorbars represent uncertainty of the fit ( fT3l) . 
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7 Summary and Discussion 



In this paper we investigated the mechanism of tunnehng via unstable semiclassical solutions 
(sphaleron-driven tunneling) which governs the processes of multidimensional tunneling at 
energies higher than some critical value Ec- There were two aspects in our study. First, we 
analyzed the experimental signatures of sphaleron-driven tunneling. These are suppression 
of the tunneling probability by the additional power of the semiclassical parameter h and 
substantial widening of the final-state distributions as compared to the case of ordinary 
barrier tunneling. 

The second aspect of this paper was related to the development of the modified semiclas- 
sical technique (the method of e-regularization) , which is applicable in the case of sphaleron- 
driven tunneling. This method is completely general; it was derived from first principles using 
the formal operations with the path integral. Similar modified technique has been imple- 
mented in several quantum mechanical [261 [IS] cind field theoretical [32] tunneling problems. 
Using the modified method, we obtained expressions for the inclusive and exclusive tunneling 
probabilities in the case of sphaleron-driven mechanism, investigated the "phase transition" 
between the cases of potential and sphaleron-driven tunneling. We also derived relation 
between the probabilities of tunneling from the low-lying and highly excited in-states. 

Our results for the power-law dependences of the semiclassical prefactor are summarized 
in Tabled! 





Potential 


Sphaleron-driven 


Inclusive, 


Ey 


~ 1 




h 


Exclusive, 


Ey 


~ 1 


h 




Inclusive, 


Ey 


~ h 






Exclusive, 


Ey 


~ h 




^3/2 



Table 1: Summary: the power-law dependences of the tunneling probability in two dimen- 
sions. 

Let us comment on the relation between the mechanism of sphaleron-driven tunneling 
and Wilkinson formula for the tunnel energy splitting [6l [Tj IHl IS]. The latter formula is 
applicable in the cases of near-integrable (as opposed to completely integrable) systems with 
double-well potentials. It is based on the following property of near-integrable dynamics: 
tunneling trajectories stemming from the wells of near-integrable system do not end up in 
opposite wells (as in the integrable case), but rather get attracted to a certain unstable 
periodic orbit0- Due to this feature the sphtting in Wilkinson formula is suppressed by the 
additional factor ^^/^ as compared to the case of completely integrable system. One observes 
that, technically, the reason for this factor is similar to that in the mechanism of sphaleron- 
driven tunneling considered in this paper. However, the two cases are physically different: 
transition to the sphaleron-driven regime is unrelated to the transition from integrable to 

^^This orbit is the intersection of the Lagrange manifolds associated with the two wells [B]. 
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non-integrable dynamics. In addition, the relevant periodic orbit in the Wilkinson formula 
is complex while the sphaleron orbit is real. 

We finish this paper with remarks on the recent observation [29j that the new tunneling 
mechanism generically leads to anomalously large times of tunneling. Indeed, the semiclas- 
sical trajectories describing sphaleron-driven transitions spend infinite time interval in the 
vicinity of the sphaleron orbit. Clearly, the time scale At of such transitions should be large, 
in particular, one expects At +00 as h —>■ 0. A rough estimate of At can be obtained as 
follows. Due to quantum fiuctuations the system cannot approach the sphaleron orbit in the 
phase space closer than at the distance determined by the uncertainty principle, ApAa; ~ h. 
The semiclassical trajectories starting in the vicinity of unstable sphaleron go away from it 
exponentially with time; thus it takes them the time At ~ log | Ap|, log | Aa;| to leave the 
sphaleron neighborhood. This translates into the characteristic lifetime of the sphaleron 
At cx I log h\ , which sets the characteristic time scale for sphaleron-driven tunneling. The 
dependence of tunneling time on h provides another possible experimental signature of the 
new tunneling mechanism. Yet more signatures can be found by analyzing the probability 
distribution over tunneling time. The modified semiclassical method proposed in this paper 
allows comprehensive study of these issues which will be published elsewhere [¥T] . 

Acknowledgments. We are indebted to F.L. Bezrukov, S.V. Demidov, D.S. Gorbunov, 
M.V. Libanov, N.S. Manton, V.V. Nesvizhevsky and V.A. Rubakov for useful and stimulating 
discussions. This work was supported in part by the RFBR grant 08-02-00768-a, Grants of 
the President of Russian Federation NS-1616.2008.2 and MK-1712.2008.2 (D.L.), Grant of the 
Russian Science Support Foundation (A.P.), the Fellowships of the "Dynasty" Foundation 
(awarded by the Scientific board of ICPFM) (D.L. and A. P.) and the Tomalla Foundation 
(S.S.). The numerical calculations were performed on the Computational cluster of the 
Theoretical division of INR RAS. 



A Semiclassical tunneling probability 

In this appendix we give details of the standard method of complex trajectories. The main 
idea of the method is presented in Sec. 14. 1[ 

Our starting point is the path integral representation (j4]) for the out-state wave function 
"^f. This representation contains two main ingredients, the in-state and the quantum 
propagator written as a path integral. In accordance with the discussion in the main body 
of the paper, the in-state has definite values of the total energy E and ^/-oscillator energy Ey. 
One writes \E'i(x, y) as a product ipx{x) ■ ipy{y), where ip^ is a plane wave with momentum 
Px^i = \/l{E — Ey) and unit fiux normalization, while '\l)y represents the semiclassical wave 
function of the oscillator with energy Ey. Combining i\}x and V'y, one obtains. 
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In this formula Py,i{y) = ^/'2Ey — uP-'xp is the y component of the momentum in the in- region 
X — oo, while 

/y 
VyAv') dy' (50) 
^2Ey/ui 

stands for the classical action in this region. Note that in Eq. fj49l) we keep only one of the 
two exponents entering the standard expression for the oscillator wave function. The reason 
is that ^j(a;, y) will be used deep inside the classically forbidden region, where the omitted 
exponent is negligiblj^. 

At small h the path integral for the quantum propagator is evaluated by the saddle-point 
technique. The result is given by the Van Vleck formula 



[dx] 



f „i5[cc(^)]/ft 
gi5[cc]/n ^ ^ 



det 



dxjdx 



1/2 



(51) 



We refer the interested reader to Ref. [H] for derivation. The formula flSTj) is written in 
terms of the semiclassical trajectory x^^\t)^ which has the meaning of a saddle-point path 
saturating the path integral. This trajectory satisfies the classical equations of motion; it 
starts from x = xi aX t = tj and arrives to x = Xf aX t = tf. Below we omit the superscript 
(s) of the semiclassical trajectory. 

We substitute the semiclassical expressions ( H9|) and (ISTl) into Eq. (jl]) and take the saddle- 
point integral over Xi. The result for the out-state wave function has the exponential form ([6]), 
where Z)~^/^ collects all prefactors including the determinant due to the saddle-point inte- 
gration; we will evaluate D below. Note that integration over Xi changes initial conditions 
for the trajectory x{t). Namely, the extremum of the leading exponent S + Bi with respect 
to Xi is achieved when 

Xi = Px,i , 2/i = PyAVi) ■ (52) 

One finds that these conditions are equivalent to the fixation of the in-state quantum num- 
bers, Eqs. dS]). 

A remark is in order. We consider the case of classically forbidden transitions which 
implies that there is no real solutions starting in the in-region with fixed E, By and arriving 
into the out-region at t = tf. Accordingly, the saddle-point trajectory x{t) is complex. 

Given the final state wave function, one evaluates the inclusive probability of tunneling 
performing the saddle-point integration over yj in Eq. ([7]). One obtains the familiar semi- 
classical formula ([T]) for the inclusive tunneling probability, where the leading exponent Bp^ 
is given by the value of the action functional (Q evaluated on the complex trajectory x{t). 
The prefactor will be discussed shortly. 

Let us comment on the final boundary conditions ([8]) obtained after integration over yf. 
One finds that all of them have different origin. Namely, the final value oi Xf is already fixed 
in the probability formula ([7]); the condition ijf = yj corresponds to the extremum of the 
leading semiclassical exponent with respect to yf. The third condition, namely, reality of 
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We assume appropriate choice of the branch ofpy^i{y), see e.g. Ref. [42] 
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i/f, follows from uniqueness of complex trajectory, which is assumecO. One also notes that 
the trajectory x(t) is real in the out-region. Indeed, Xf, yj, ijf are real due to the boundary 
conditions at t = tj, while e M due to conservation of total real energy E. 

We finish the discussion of the complex trajectory x{t) by remarking shortly on the 
important issue of non-trivial contour in complex time [22l |38]- We noted already that x{t) 
is real in the out-region. Solving the classical equations of motion backwards in time, one 
concludes that at t G M the trajectory is real as well. Thus, x{t) in real time corresponds 
to classically allowed reflection from the barrier; clearly this solution is not relevant for the 
description of tunneling. One observes, however, that the semiclassical trajectory has a 
branch point in the complex time plane, see Fig. [TH The solution describing tunneling is 
obtained along the contour winding around this branch point (the contour A'ABCD in the 
figure) . 




Re t 



Figure 14: The tunneling solution is obtained along the contour A'ABCD in the complex 
time plane. The branch point of the solution is shown by the cross. 

Let us evaluate the prefactor Apot- There are two non-trivial contributions into Apot 
coming from the saddle-point integrations in Eqs. (jlj) and (JTj). We start with the prefactor 

I. One finds that 



D Qf ^Jie out-state wave function, Eq. (p 

dXf 



D 



2tt 2 
— i/iXi X {2mn) 

UJ 



det 



dxi 



-1 



■det 



d^{S + Bi 



dx. 



X f =const 



271 . . . d\S + Bi 
— yiXi- det — — — — 

CU OXiOXf 



(53) 



:cf =const 



where the three multipliers in the first line correspond respectively to the prefactors of the 
in-state (H9|) . propagator (15T]) and the determinant due to the saddle-point integration over 
initial coordinates Xi. In the second equality we combined the multipliers. 

One substitutes the explicit form of Bi into Eq. (153|1 and represents the factor yi Xi as 
the determinant of diagonal matrix. The result is 



2vr , 
L) = _ det 

UJ 



dxj 
dx 



f 



\ dxi 



dx 



(54) 



^^The condition yf — should be relaxed if several complex trajectories contribute into the out-state of 
the process. In this case the trajectories with complex yf give rise to interference terms in the tunneling 
probability. 
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This last determinant can be evaluated by considering the set of linear perturbations Sx{t), 
which satisfy the linearized equation of motion fllOl) in the background of complex trajectory 
x{t). The two basic perturbations will be particularly important, il^^^\t) = dx{t)/dxf and 
il)^'^\t) = dx{t)/dyf, where the derivatives are taken a.t Xf = const. One can explicitly check 
that satisfy Eq. (fTOl). 

Consider linear energy increment due to perturbation 6x{t), 



6E[6x] = Xi6xi + 6Ey[6x] 



(55) 



where the increment SEy[Sx] of the initial oscillator energy is given by Eq. ( |T2l) . Clearly, 6E 
is conserved; in particular, it can be computed at t = tf. Using the energy increments, one 
rewrites the determinant dSH) as 



2vr , 
D = —det 



2ti 
27r 



det 



bEy\^^'\ 

2// 

(2) 



(56) 



where in the last two equalities we added the second row to the first, computed explic- 
itly bE\^^^\ =3;/, bE\ilP'\ = ijf and used linearity of 5Ey. Let us introduce the linear 
combination 

5x^^\t) = -yf/xf^P^^\t)+iP^^\t) . 
This vector satisfies the linearized equations of motion with the Cauchy data at t = t/. 



Note that SE[Sx 



(1)1 



0. Using 6x^^\ one finds. 



D 



271 Xf 



UJ 



6Ey[6x 



5xf 



(1)1 



. 



(57) 



(58) 



So, we have obtained the final state prefactor D^^^"^ in Eq. ([6]). 

Let us proceed with the probability prefactor Apot- Taking the saddle-point integral in 
Eq. ([7]), one arrives at the following formula. 



A 



pot 



Xf\fnh 



\DyD' 

where the factor coming from the integration is 

.9yf 



(59) 



D' = Im 



(60) 



E, Ey , X j^=const 
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One evaluates D' considering the perturbation 



Xf,Xf=const 



It satisfies the hnearized equations of motion fllOl) with the final Cauchy data 

6xf = i-Lu'yf/xf, 0) , 6xf = (0, 1) , (61) 

so that = 0. 

One notices that 6x^^\t) and Sx^'^^t) do not change the values of and E. Moreover, 
these two are the only linearly independent perturbations which have this property. Thus, 
the perturbation 

dx{t) 



Pit) 



(62) 



E, Ey, xf=const 

is their linear combination, p(t) = aSx^^\t) + pSx^'^\t). One has, 

Pyitf) = l ^ /? = 1, 

5Ey[p]=a6Ey[5x^'^]+P6Ey[5x^^^]=0 ^ = -^^j^ • 
From Eqs. (I60l) . (|62|) we deduce the formula 

D = Imp,(t,) = Ima = prj^^(T)]p • (63) 

Substituting expressions (158|) and (l63l) into Eq. (l59l) . one obtains Eq. (fTTi) . 

We finally rewrite the Cauchy data (1571) . (!6T|) for perturbations in canonically covariant 
form. The new conditions are: 

(i) perturbations dx^"^^ and their momenta ^i*-"-* are real at t = tj; 

(ii) they do not perturb the total energy, 

(iii) their norm is fixed by ^l{5x^-^\ 5x^'^^) = 1, where = dpx A dx + dpy A dy is the 
symplectic form. 

Note that the conditions (i) — (iii) do not completely fix the perturbations 5x^^\ Namely, 
a linear admixture of the perturbation xi^) = x(t), if added to Sx^"'\ does not disturb (i) — 
(iii) and the value of the pref actor Apot- One shows this exploiting the properties of Xy 

6E[x] = SEyix] = , n{6x^'\x) = ^{Sx^'\x) = . 

In practice one fixes the x-degeneracy by supplying some additional Cauchy datum. 

It is straightforward to check that the conditions (i) — (iii) follow from Eqs. (157|) . (!6T!) . In 
fact, they are also sufficient^ i.e. any pair of perturbations Sx^^"* and 5cc^^-' satisfying (i) — (iii) 
can be used in the prefactor formula (fTTj) . One shows this explicitly by decomposing the new 
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perturbations 5x^'^^ in the basis of the old perturbations, substituting them into Eq. (fTT!) 
and using the properties (i) — (iii). 

Let us note the following property of the formula fill I) for Apot- The perturbations 6x^'"'\t) 
leave unchanged all boundary conditions of the complex trajectory except for one, fixation 
of Ey. It is precisely the change in this datum which enters the formula for the prefactor. In 
the limit of separable system, when Ey becomes a conserved quantity, 5Ey[5x^'^^] becomes 
real due to the condition (i) and Ap^ tends to infinity. This indicates the change of the 
/i-dependence of the prefactor in this limit, cf. Refs. P,[H1[S]- 

B Numerical method for semiclassical calculations 

Here we describe the numerical method of finding the trajectories x{t). The method has 
two useful properties. First, it is applicable to systems with arbitrary number of degrees 
of freedom, up to the field theory case {Af = oo). Second, it naturally incorporates the 
modified semiclassical technique of Sec. 14.21 Originally, the method was proposed in the field 
theoretical context ^51 [32] . It was adapted for quantum mechanical problems in Ref . [38] . 

We compute x{t) by solving numerically the classical equations of motion with the bound- 
ary conditions (ISl) and (IE]) imposed at t = tj and t = tj respectively. To this end we introduce 
non-uniform lattice {tfc, k = 1 . . .Nk}-, where ti = ti, t^^ = tf, and discretize the classical 
equations of motion and boundary conditions in a straightforward manner. We find that 
the second-order discretization works well enough. As mentioned in appendix [Aj the time 
variable t runs along the contour in complex time plane, see Fig. [TH Accordingly, the sites 
tk of the lattice belong to this contour. 

After discretization one obtains the system o{2 x Nj. complex non-linear algebraic equa- 
tions for the same number of unknowns Xj. = x{tk). Let us denote the unknowns collectively 
by Za, where a = 1 . . . 2Nk, and equations by J-'b{z) = 0. We solve the equations by the 
Newton-Raphson iterative method, see e.g. Ref. |16|. In this method one starts with some 
approximation z = z^^^ for the solution. Then the approximation is repeatedly refined by 
finding corrections Az from the system of linear equations 

+ Az) ^ V ^(zW) Az, + :fM'^) = . 

Note that the coefficient matrix of this system is block three-diagonal and can be inverted 
efficiently. At the end of each iteration one redefines the approximation, z^^^ — ^ z^^^ + Az. 
After 3 — 6 iterations the method converges provided the original approximation was good 
enough. 

The drawback of the Newton-Raphson method is the small radius of convergence. Namely, 
the method does not produce correct solution unless the approximation z^^^ is sufficiently 
close to it. We solve this difficulty in the following way. Suppose the solution is known 
for some values {E, Ey) of the in-state quantum numbers. Then, one finds solution at 
{E+AE, Ey + AEy) by the Newton-Raphson iterations, with the original solution at {E, Ey) 
serving as the initial approximation. Using this approach one walks in the {E, Ey) plane by 
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changing the values of quantum numbers in small steps and finding the respective solutions. 
Moreover, in this method one can gradually change any parameter of the problem including 
the regularization parameter e of Sec. 14. 2[ 

The final ingredient of our numerical procedure is the method of finding the semiclassical 
solution at some special values of E and Ey. We start the procedure by computing the 
instanton trajectory which describes tunneling at E = Ey = 0. This trajectory is real in 
Euclidean time and can be obtained by minimization of Euclidean actioio- After finding 
the Euclidean instanton, one bends the time contour in a way shown in Fig. [T3] and finds 
the in- and out- parts of the trajectory by solving the Cauchy problem from points B and C 
of the contour. 

Note that there are several methods of starting the numerical procedure. In particular, 
a greater class of Euclidean solutions (periodic instantons) can be used for that purpose, see 
Ref. [38] • In the cases when neither instanton nor periodic instantons exist one can exploit 
classical over-barrier solutions [iBj. 



C Saddle— point integrals in the modified method 

The integral over e in Eq. (fT8l) is evaluated as follows. One notes that by construction the 
trajectory x^{t) extremizes the functional + Bi. Thus, 

^(S.|..l + = ^(%±^ . ^ + ^(%±^ = ir.„, w , (64) 

de ox^ de oe 

where in the second equahty we used 5{S^ + B^jbx^ = 0. Using this relation, one finds that 
the saddle point of the leading exponent in Eq. (ITSl) with respect to e is achieved when 

TrntiXe] = r . (65) 

The result for the final wave function is 



Jo y/2nhD^ V dr 

Note that the saddle-point value of e does not need to be purely imaginary. 

The probability formula ([7j) involves, besides the complex conjugate out-state ^j. 
One derives the analog of Eq. fl66l) for \I'j considering the path integral, which is complex 
conjugate to Eq. (jl]). We substitute expressions for j and \I'j into Eq. ([7j) and obtain, 

V= f dVf J . ^iiS4^.]+B,[a^.]-ieT-S_A'^_,,]-Bda^_^,]-te'T')/h /g-,N 

J Jo 2nhVD,D_,, V drdr' ' ^ ^ 

where the integral over r' comes from the conjugate out-state. Note the opposite signs of e 
and e' in Eq. fl67|) : the difference is related to the fact that the modified action f|T7|) depends 



•"Say, with the algorithm of conjugate gradients [46| . 
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on the combination ie which changes the sign under complex conjugation. As a consequence, 
the saddle-point condition for e' reads Tint[x_^'] = r', cf. Eq. (^5^. 

The integral over yf in Eq. (I^7|) is evaluated in the same way as in appendix |X1 Below 
we consider the integrals with respect to the interaction times r and r'. One changes the 
integration variables to r_ = r — r' and r+ = (r + t')/2. We noted in the main body of the 
paper that fixing r+ one stabilizes both trajectories and x^/. The integral over r_ is taken 
by the saddle-point method. One finds the extremum of the leading exponent with respect 
to r_, 

e' = e , (68) 

where the relation flM|) was used. Note that after the integration over r_ the value of e is 
defined by the implicit relation 

Tint[x,]+T,r,t[X-,]=2T+ . (69) 

One finds that for real r+ the solution e = e(r+) of this equation is real. To show this we 
assume that the complex trajectory x^{t) is unique. Then, it is straightforward to check that 
the semiclassical equations following from implj^that x* = x_^*. Therefore the l.h.s. of 
Eq. (169|) is real for real e, and so is the function T+(e). This entails the reality of the inverse 
function e(r+). The condition e = e* and Eq. (!69|) are equivalent to Eqs. (fT9|) of Sec. 14.21 
The result of integration over r_ is given in the main body of the paper, Eq. (120|) . 



D Evolution near the sphaleron 

In this appendix we study the evolution of the system in the vicinity of the sphaleron orbit. 
The aim of this analysis is to extract the behavior of the suppression exponent and prefactor 
in the regime when the tunneling trajectory spends a long time near the sphaleron. This 
is the case for the (regularized) inclusive trajectories at e -C 1 (Sec. 14.21) and exclusive 
trajectories with large topological numbers m (Sec. E]). 

Let us start with the limit e +0 in the modified expressions 0221) . We work in the 
approximation of small sphaleron amplitude. Though this approximation is justified only at 
energies slightly exceeding the minimum height Vq of the potential barrier, we believe that 
the qualitative features of and Apot,e remain the same at higher energies. 

In small vicinity of the saddle point the potential is approximated by 

V{x, y) = Vo + {ujIxI - Lj^_xl)/2 + 0(x=^) , 

where the Cartesian coordinates x+ and X- run along the stable and unstable directions of 
the potential, while uj± represent the respective frequencies!^ The sphaleron orbit describes 

^^If the trajectory is not unique, the relation a;* = aj-^. is no longer valid for the terms in the tunneling 
probability which account for the interference between different trajectories. We do not consider interference 
effects in the present paper. 

^^In the model (H]) Vb = 1, lu± = ±{lu^/2 — 1) + ■\/a;'*/4 + f and the coordinates {x+, xJ) are rotated with 
respect to (x, y) by the angle a = iarcctg(ijj^/2). 
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periodic oscillations along 



where a+ is related to the sphaleron energy. 

Now consider the modified complex trajectory x^{t). At small e it has two distinctive 
parts corresponding to the two stages of the tunneling process. First, the trajectory arrives 
into the vicinity of the sphaleron orbit. Second, it leaves the sphaleron and evolves into the 
out-region. In the vicinity of the saddle point one writes 

= x'f{t) , = a_e-"-* + ea_e+"-* . (70) 

In writing down Eq. fITOl) we took into account two facts. First, at e = the exponentially 
growing term in the equation for x^^^ vanishes and the trajectory stays forever in the vicinity 
of the sphaleron. Second, this term is proportional to e due to the linear dependence of the 
modified equations of motion on e. At small e Eq. flTUl) describes long intermediate stage of 
modified evolution near the sphaleron. 

The exponentially growing term in Eq. (!70l) destroys the sphaleron orbit within the time 
interval 

r+ = -— lne + 0(l) . (71) 
Using the Legendre transformation (pT]) we find, 

= F,ph + / dT+ 2e(r+) = F^ph + — + 0{e^) . 



Therefore, tends linearly to its limiting value. 

Consider now the limit e +0 of the prefactor. One finds Apot,e by considering linear 
perturbations 5x'^'^\t) in the background of the modified trajectory. Namely, one fixes the 
Cauchy data for Sx'^'^^ at t = tj and evolves them backwards in time. As one reaches the 
stage of near-sphaleron evolution, the perturbations start growing: they contain the part 
ri{t) ~ e""^-* which grows exponentially as t decreases. One writes, 

5aj(") (t) = ■ rt{t) + 6xlll{t) , (72) 

where the last term stay bounded as t decreases. The coefficients d^"^ are real due to the final 
Cauchy data. Consider now the perturbations 5x^"'\t) at t = ti. One observes that the two 
terms in Eq. ( 1721) behave differently in the limit e +0. Namely, r]{ti) ~ e^-'^+ ~ 0(l/e), 
while the second term is finite, 6xi^l{ti) ~ 0(1). Using this dependence, on finds from 
Eq. (HH) that Apot,e ~ 0{e^/^). This fact and Eq. (^^ imply that the limit fl22bl) for Asph 
exists. 

Our next goal is to prove Eq. fl35l) . In what follows we drop the assumption of small 
sphaleron amplitude. Consider the sphaleron orbit Xsph{t)- This orbit is a periodic solution 
of the equations of motion with period Tgph, Xsph{t+Tsph) = Xsph{t)- It is completely specified 
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by two parameters, the total energy E and time origin to. A small perturbation Sx around 
Xgph satisfies the equation 

Sit + V"{xsph)5x = . (73) 

This equation has two obvious solutions, 6xi = dXgph/dE, 6x2 = Xgph', these are the deriva- 
tives of the sphaleron orbit with respect to its two parameters. The two remaining solutions 
of Eq. (1751) describe the formation and decay of the sphaleron. According to Floquet theorem, 

Sx.{t) = Sx^{t)e-^' , Sx^t) = 5x+{t)e^' , (74) 

where Sx^{t), 6x^{t) are periodic functions with period Tgph and f3 is the Lyapunov exponent. 
In what follows we restrict our attention to the perturbations ([71]): we omit the mode 6x1 
because we are interested in perturbations preserving the total energy; the mode 6x2 is 
removed by the trivial time shift. 

The tunneling trajectory x{t) with energy E has the following form in the vicinity of the 
sphaleron, 

x{t) 

'^sph 

(t) + 6x^{t) +C+6x+{t) . 

The coefficients C_, C+ completely parameterize the trajectory. Varying these coefficients 
one goes over the possible values of the initial and final oscillator energies, Ey and i?^. Let 
us consider the suppression F calculated on the tunneling trajectory as a function of Ey, 
C+. The choice (7+ = corresponds to the trajectory which stays at the sphaleron forever. 
Clearly, F{Ey, C+ = 0) = Fgph- At small but non-zero values of C+ one has, 

+ 0{Cl) . (75) 

Ey=const 

Consider now two exclusive trajectories which have the same E^ and the topological num- 
ber^m and m + 2. As discussed in Sec. |5l the (m + 2)-th trajectory performs one additional 
oscillation in the vicinity of the sphaleron orbit as compared to the m-th trajectory. This 
implies that the coefficients and corresponding to these trajectories are related 

by (see Eq. dUD) 

= e-'^^-'^Cf ) . (76) 

Note that the Lyapunov exponent (3 in this expression depends on the properties of the 
sphaleron orbit only; in particular, it is independent of the final oscillator energy. Substitut- 
ing the relation fl76l) into Eq. fl75|) yields Eq. fl35|) . 

Using the results of this appendix, one estimates the parameter j3 in Eq. fl5^ : (3 = 
TsphP ~ 27rcj_/cj+. In the model f3 ^ 2A. 

E Implementation of the uniform formula 

The quantities entering the uniform correction factors fl28|) are computed as follows. 
^■^We restrict our attention to the trajectories of the first class in the terminology of Sec. [5l 



^sph + 



dF 



9C, 
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Consider first Aipot- One obtains the saddle-point value Wg by taking the integral in 
Eq. fl24p where it is convenient to change the integration variable to e, 



Ws = -i= [ de' J-dr'^/de' ■ Apot,e' ■ (77) 

Here e = ej < corresponds to r+ = +00 and we used the fact that the saddle point is 
achieved at e = 0. Note that the integrand /(e) = —drj^/de ■ Apot,t of Eq. ( 1771) is singular 
at e = ej, since it contains the derivative of r+ in the nominator. We plot this integrand 
(points in Fig. [T5k) for several values of energy and Ey = 0.05. We find that the function 
/(e) is well fitted by the formula /(e) ^ A/{e — e^)^/^ + B (lines in Fig. [T5b). We exploited 
this fact in the numerical computation of the integral in Eq. (1771) . The relative numerical 
error of Ws was always smaller than 10^^. The result for the saddle point w = Ws{E) is 
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Figure 15: (a) The integrand in Eq. fl77j) plotted as function of e for several values of total 
energy E < Ec{Ey) and Ey = 0.05. The limiting values e = e^ are shown by the vertical 
dotted lines, (b) The saddle point Ws{E) at Ey = 0.05. 

plotted in Fig. [T5b . As expected, Ws is positive, decreases with energy and reaches Wg = 
at E = E^{Ey) ^ 1.096. 

The second derivative of the suppression exponent entering Eq. (I28ap is computed using 
the formula 

F"{w.) = ^ , (78) 

which follows from the definition of w, Eq. (12^ . Using Wg and F"{ws), one finds the argument 
Xpot of the Fresnel integral and thus J^pot, see Eq. fl28ap . 

Consider now Aigph- From Eqs. (1241) . ( |2T]) one derives in a straightforward way the 
following expressions for the derivatives of the suppression exponent at w = 0, 



F'{0) 



h 



A 



sph 



h 



A 



sph,€ 



(79) 



e=0 
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Here we denote by Agph^^ the value of the r.h.s. of Eq. fl22b|) at finite e. Combining the 
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Figure 16: Suppression exponent as function of \ogAsph,e- The graph is plotted at i^^ = 
1.12 > E^{Ey), Ey = 0.05. 

multipliers fl79p . one obtains for the variable entering Eq. (]28bl) . 



Note that x^p/j is real and positive. One finds it by plotting F^ as a function of log Agph^e at 
small e (Fig. [TBI) and fitting the graph with the linear function. 

F Semiclassical probability at small Ey 

We showed in Sec. [6] that the probability of tunneling from the ground state of ?/-oscillator 
is given (up to the overall factor (vr/e)^/^) by the standard semiclassical expression, where 
one should use Ey = huj/2. Let us determine explicitly the ^-dependence of this expression, 
i.e. extract the leading exponent and prefactor of the probability. 

First, we rewrite Fpot using the parameters T and 6. Namely, we evaluate the in-state 
term Bi by taking the integral in Eq. (l50l) . 



where in the last equality the asymptotic form (l45l) . (I46l) of the trajectory was used. For the 
suppression exponent ([9]) one obtains. 



sph 



1 dF, 

2h dlogAsph,e 



(80) 



2Imi?j[a;] = Im {yii/i + XiXi) + Im \XiXi + i2Ey/uj) arccos{iuyi/ ^JlE.y) 
= lm{yiyi + XiXi) - 2ET - Ey9/uj , 



Fpot = 21mS[x] - 2ET - Ey9/u , 



(81) 



where 




(82) 
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is the classical action integrated by parts. 

The second observation is as follows. Consider the differential of the action fl82l) with 
respect to the in-state quantum numbers E, Ey. One writes, 

d2\TsiS = dim {2S + XiXi + yiiji) = Im {xidxi - Xidxi + yidiji - ijidyi) 
= 2EdT + Eyde/u . 

This equation together with Eq. ( IHTl) show that the parameters (T, 9) are related to {E, Ey) 
by the Legendre transformation. Consequently, 

dEpot = -2TdE - OdEy/uj . 

Thus, (— 2T) and {—O/uj) are equal to the derivatives of the suppression exponent with 
respect to E and Ey. 

Using the above observation, we expand the suppression exponent around the point 

Ey = Q, 

1 



Ppot Ppot 



9{E'y)dE'y. 



Ey=0 LU 

The function 9{Ey) at small Ey is determined by noting that 

Ey = 2uj^\a\'^e-'^'^^-^ , 6 = -\n{2Ey/uj) + 60 + 0{Ey) , (83) 

where we used the fact that \a\ has a well-defined limit as Ey ^ 0. This yields the expression 
(HSj) for the suppression exponent. Finally, substituting Eq. (HSj) into Eq. ([T]) and recalling 
the additional factor (vr/e)^/^ one arrives at the expressions fHTl) . 
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